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Abstract 

The  long  wavelength  modes  of  excitation  of  a  two 
component  plasma  In  a  steady  magnetic  field  are  examined. 
Two  linearized  Boltzmann  equations  are  given  with  collision 
terms  whoch  are  coupled  through  the  difference  in  tempera- 
tures and  difference  in  velocities  of  the  two  gases.   A 
formal  means  of  classification  of  phenomena  is  described  in 
terms  of  the  nature  of  the  roots  about   k  =  0.   Two  types 
of  behavior  are  uncovered:   MHD,  Magnet ohydrodynamic s ,  which 
includes  finite  phase  speed  phenomena  near  k  =  0;  and  PEM, 
Plasma  Electromagnetics,  which  includes  infinite  phase  speed 
phenomena  near  k  =  0.   The  disslpative  effects  of  collisions 
are  included.   In  the  limit  of  vanishing  collision  frequency, 
roots  previously  obtained  are  recaptured.   The  relevance  of 
the  pertinent  domains  are  discussed  and  the  complex  inter- 
play between  the  "fast"  and  "Alfven"  modes  of  MHD  and  the 
plasma -magnetic  modes  of  PEM  are  demonstrated.   Equations 
macroscopic  in  appearance  are  derived  which  include  the 
effects  of  the  initial  configuration.   In  the  limit  of  large 
collision  frequency  these  equations  reduce  to  standard  forms. 
The  dispersion  of  the  N   order  Larmor  resonance  is  given 
which  includes  the  effect  of  the  mass  ratio.   A  discussion 
of  the  transfer  equations  of  a  plasma  is  included. 
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Errata 

1)  p.  25    eq.  23:   nxT  should  read  hkT  . 

2)  p.  44    eq.  14:   exponential  should  read 
exp  .^  aa '  -  — - —  [  sin  a  -  sln(a  +  a '  ] 


3)  p.  45  eq.  20:   2^^^  should  read  e^^^  . 

4)  p.  58  line  9:      Uj_      should  read  Uj^  . 

5)  p.  70  line  3:   3a^+  1  should  read  3a^-  1   . 

6)  p.  73  add   -2ik^  x  B[S^]^°^(vn  +  ■^)      to  eq.  34  . 

7)  p.  81  eq.  7:   M-^J  =  e/c^   should  read  \x^J  +   E/c^   . 

8)  p,  91  eq.  28:   l/lo)  should  read  l/icn  . 

9)  p.  110  footnote  page  reference  is  p.  38 

line  16:   change  p   to  ^   . 

10)  p,  130  In  the  equation  for  I»  change  2in  sin  ■^  to 

2iy  sin  1/    . 

11)  p.  5  General  Notation:   X  =  X-X  should  read  AX  =  X-X. 

12)  p.  88  eq.  19:   first  term  should  read   -icokP^U^°'   . 

13)  p.  131  omit  third  equation. 


(■  y 


:i    .q      (Cl 


MP-22  N^°-9^^2 

PROPAGATION  OF  LONG  WAVELENGTH  DISTURBANCES  IN  A  PLASMA 
R.  L.  Llboff 

Errata 


1)  P- 


5       General  Notation:   X  =  X-X  should  read  AX  =  X-X. 


2)  p.  14      Diagram:   Lower  left  corner;  change  cd  =  0  to 

(P0)(2^ 

3)  p.  17      Line  13:   change  to:  "...,  then  these  roots  collapse 

2     2  „ 
to  the  four  plasma  oscillations  co^  =  a>^- 


k)      Sec.  B-I   eq.  10 
eq.  19 


change  the  first  n  to  n. 
change  W  to  W. 
eq.  23:   change  nKT  to  nKT. 

5)  Sec.B-II   eq.  2:   change  H^  ^  to  H^   .  ^ 

eq  10:   The  right  hand  side  includes  the  factor 

eq.  23:   change  fi  to  Q,. 

6)  Sec.C-I    eq.  3:   In  equation  for  h-,^,  change  u  to  u. 

eq.  1^:   A  minus  sign  should  follow  aa •   . 

eq.  20:   change  2  ^  to  e 

eq.  32  and  3^:   change  4^^  to  c^^  . 

7)  Sec.  C -II   eq.  1:   change  h  to  -h'. 

eq.  10:   change  ^'B  fi^     to  4-Bxtr^. 

eq.  18:   change  |-  factor  to  1. 

p.  58,  line  9:   change  u^  =  F  Au  to  u_^  =  F  Au. 

eq.  30:   should  read:   -  =  ^^^T  +  r""  +  =~   * 


I 
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I)   Sec.C-III 


eq. 


N 
Change  cos   to  cos 


.M 


Appendix  C. 
Appendix  D. 


eq. 

27.1: 

eq. 

30: 

eq. 

34: 

eq. 

39: 

9) 

Sec. C -IV 

eq. 

7:      c 

eq. 

19: 

eq. 

28: 

eq. 

30: 

eq. 

38: 

eqs.  13:   should  include  i^^    factor  in  integrands. 

2  2  2 

eq.  25:   First  k  should  be  k  and  the  second  k 

2        " 
should  be  k^^  . 

■^  should  read  -^r  . 
oT  ot 

2      -^2 

change  C   to  C  . 

add  the  term  -2ik^  x.  B[S^] ^°^ (vn  +  "/) . 
Insert  a  plus  sign  between  (vn  +  '7^)   and  h 


eq.  7:   change  second  equals  sign  to  a  plus  sign. 

first  term  should  read   -icok  .J-^^U^^' 

change   1  to  i. 

change     J><;B     to     JX,B. 

change     Q     to     Q. 

eqs.  40:   Left  sides  should  include  the  factor  2 

and  the  radicals  should  be  multiplied  by  i 

2  2 
p.  103:   line  seven,  omit  k  A  . 

eq.  64.1:   all  speeds  should  be  squared. 

eq.  69:   change  fi  to  f2  . 

p.  13O:   In  the  equation  for  I„   change  2in  to  21^ 
p.  131 :   Omit  third  equation. 
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Propagation  of  Long  Wavelength  Disturbances  In  a  Plasma 


Introduction 


The  problem  of  uncovering  the  modes  of  excitation  of  a 
plasma  In  a  uniform  magnetic  field  has  been  extensively  studied 
In  the  literature.   For  the  most  part  these  endeavors  were 
constructed  about  a  macroscopic,  phenomenologlcal  theory.   In 
addition,  one  Is  able  to  find  many  problems  which  yield  solu- 
tions under  special  constraints,  such  as  on  the  direction  of 
propagation  of  the  disturbances;  statements  pertaining  to  the 
vanishing  of  certain  macroscopic  variables  (e.g.,  low  tempera- 
ture approximation,  charge  neutral  assumption,  etc.)  and  ad  hoc 
assumptions  relating  the  higher  moments  of  the  distribution 
function  to  Its  lower  moments. 

These  Investigations  lead  to  distinct  categories,  e.g., 
magnetohydrodynamlcs,  '    '   cooperative  phenomena   (plasma 

oscillations),  magnetic  resonance,  '''   cold  plasma  theory,^' 

11  12 
colllslonless  limit  (Vlasov  equation),   '    etc. 

In  this  analysis  we  have  examined  the  general  problem 
of  uncovering  the  modes  of  development  of  a  plasma  In  a 
uniform  magnetic  field.   The  plasma  Is  displaced  from  an 
equilibrium  state  which  Is  described  by  two  absolute  Maxwelllans 
with  equal  temperatures,  equal  number  densities,  and  zero  macro- 
scopic velocities,  together  with  constant  magnetic  field  and 
zero  electric  field.   The  state  of  the  system  Is  subsequently 
described  by  the  perturbation  of  these  variables  away  from 
equilibrium. 


The  principle  concern  has  been  to  present  a  self-contained 
approach  to  the  problem,  starting  from  microscopic  origins,  and 
in  this  manner  to:   (l)  find  all  the  modes  of  excitation,  (2) 
exhibit  a  formal  ordering  of  phenomena  into  distinct  categories, 
(3)  uncover  the  connections  between  these  modes,  especially  with 
respect  to  the  collapse  of  entirely  different  domains  to  a 
common  class  of  phenomena,  and  (4)  follow  the  transition  from 
microscopic  to  the  macroscopic  description. 

We  are  primarily  interested  in  long  wavelength  disturbances j 
i.e.,  the  technique  of  solution  involves  expansions  about  small 
wave  number.   The  equations  of  motion  are  two  Boltzmann  equa- 
tions and  Maxwell's  equations  which  are  coupled  through  the 
charge  density  and  the  current.   These  equations  are  perturbed 
away  from  the  above  stated  equilibrium  state,  and  the  resulting 
system  is  linearized  in  the  perturbation  quantities.   The 
linearized  collision  integrals  are  replaced  by  forms  which  yield 
the  standard  macroscopic  equations,  and  are  consistent  with 
basic  relaxation  arguments.   In  addition  these  forms  are  shown 
to  follow  quite  naturally  from  an  expansion  of  the  linearized 
collision  integrals  in  terms  of  their  own  eigenfunctions,  or 
through  an  itteration  procedure  of  the  full  Boltzmann  equations 
about  local  Maxwellian  states,  (cf.  Appendix  A).   To  uncover 
the  subsequent  development  of  the  system  a  space-time  Fourier- 
Laplace  transform  is  introduced  which  imm.ediately  permits  the 
distribution  function  to  be  written  as  a  linear  functional  over 
the  first  five  moments  of  both  gases  (viz.,  n,n;  u,u; 
T,T).  '      '  -^   If  one  then  operates  on  the  distribution  function 
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written  In  this  manner  with  the  related  first  five  moment 
Integral  operators,  one  obtains  a  linear  homogeneous  set  of 
algebraic  equations.   These  equations,  together  with  Maxwell's 
equations,  which  relate  the  electric  field  to  the  current  and 
charge  density,  constitute  a  self-consistent  set  of  equations 
in  the  thirteen  transform  variables  which  are  the  transforms 
of  the  two  number  densities,  the  two  macroscopic  velocities, 
the  two  temperatures,  and  the  electric  field.   The  coefficient 
matrix  is  a  transcendental  function  of  the  frequency  cu  and 
wavelength  k    of  the  perturbed  state.   For  the  most  part, 
the  relations  we  find  are  asym.ptotic  to  the  exact  form.   A 
non-trivial  state  will  exist  only  for  those  values  of  k  and 
CD  which  make  this  matrix  singular.   It  follows  that  the  problem 
of  uncovering  the  modes  of  development  is  essentially  an  eigen- 
value problem.   The  related   13  x  1   column  vector  determines 
the  state  of  the  system.   For  convenience  of  reference  we  will 
call  the  said  matrix  the  master  matrix  of  this  problem  and  the 
determinant  will  similarly  be  referred  to  as  the  master  deter- 
minant , 

The  overall  analysis  is  in  two  stages.   Firstly,  the  above 
stated  formal  expansions  are  performed.   In  this  procedure  the 
parameter   ck  (c   is  the  speed  of  light)  is  treated  as  an 
Independent  variable.   Specifically,  the  lim.lt   k  — >  0   is 
taken  along  the  path   ck  =  constant.   To  terms  of  order  k 
the  master  determinant  is   l4    order  in  go.   These  roots  sepa- 
rate into  two  groups  which  are  categorized  in  terms  of  the 
behavior  of  the  i-oots  about   k  =  0.   In  the  second  stage,  the 
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inversion  equations  in  (x,t)  space  are  examined  in  order  to 
discover  the  domain  of  relevance  of  the  related  classes  of 
phenomena.   In  this  manner  we  find  that  the  PEM  (Plasma  Electro- 
magnetic) limit  is  characteristic  of  a  low  temperature  plasma, 
and,  furthermore,  that  within  this  limit  the  variation  of 
number  density  generates  an  entire  spectrum  of  phenomena,  ranging 
from  a  coupling  to  pressureless  MHD   (Magnet ohydrodynamics  plus 
displacem.ent  current)  in  the  infinitely  dense  region  to  the 
vacuum  electromagnetic  and  pure  Larmor  roots  at  the  opposite 
extreme. 

The  generalized  macroscopic  equations  which  emerge  (cf. 
Eqs.  C-IV,  1-10)  exhibit  the  following  properties:   (l)  The 
momentum  equations  contain  non-local -in-time  pressure  terms, 
which  in  the  limit  of  large  collision  frequency  reduce  to 
standard  forms.   (2)   The  source  terms  of  the  momentum  equa- 
tions include  difference   in  temperature  contributions,  in 
addition  to  difference  in  flow  terms  suggested  by  previous 
authors.  '      '  (5)   The  energy  equations  contain  non- 

isotropic  heat  conductivities  which  are  functions  of  the 
temperatures  of  both  gases.   In  the  limit  of  equal  masses  and 
temperatures  and  large  collision  frequencies  these  forms 
reduce  to  the  classical  Chapman -Cowling  result.  ^   (4)   The 
energy  equations  contain  difference  in  temperature  sources  which 
contribute  to  the  time  rate  of  change  of  the  individual  entropies, 
so  that  to  lowest  order  the  component  gases  are  non-adiabatic . 

The  classification  of  phenomena  into  two  groups  as  mentioned 
above  is  specified  according  to  the  vanishing  or  non-vanishing 
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of  the  leading  term  in  the  expansion  of  the  dispersion  relation 
about   k  =  0,   viz.,   cd  =  cd^  '   +   kco^  '   +    ,..    .   The  roots  of  the 
former  class  are  indicative  of  finite  wave  speed  phenomena  (MHD), 
while  those  of  the  later  class  (PEM)  include  infinite  phase 
speed  phenomena.   In  the  special  case  of  vanishing  collision 
frequency,  the  PEM  region  reduces  to  the  cold  plasma  domain;'  ' 
This  terminology  is  motivated  by  the  observation  that  the 
relevant  momentum  equation  follows  from  the  more  general  one  in 
the  limit  of  vanishing  temperature,  i.e.,  [kT  /MCv]  «  1,  where 
K   is  Boltzmann's  constant,   M  is  the  mass,   C   is  the  thermal 
speed,  and  v   is  the  speed  a/k.   In  the  PEM  region  v  — >  <», 
as   k  — >  0  and  the  inequality  is  naturally  established.   A 
similar  observation  is  evident  with  respect  to  the  MHD  region 
wherein  one  notices  that  the  relevant  current  equation  follows 

from  the  more  general  one  in  the  limit  of  large  number  density, 

2     2 
or,  more  precisely:   cu  »  cd  ,   where  cu   is  the  plasma 

frequency,  and  co  is  the  frequency  of  the  disturbance.   In  the 
MHD  region,   co  — >  0  as  k  — >  0,   and  this  requirment  is  auto- 
matically satisfied. 

The  dispersion  relation  of  MHD  will  be  shown  to  be  6  '  order 
in  CO,   while  the  dispersion  relation  of  PEM  will  be  seen  to  be 
10   order  in  cd.   Including  the  two  distinct  temperature  decays, 
this  complement  of  modes  includes  eighteen  roots.   Since  the 
original  master  determinant  is  l4    order  in  o),   it  follows 
that  four  of  the  six  MHD  roots  and  four  of  the  ten  PEM  roots 
are  different  forms  of  the  same  modes.   These  four  comm.on  modes 
are  shown  to  be  the  doubly  degenerate  cold  fast  waves  and  the 
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two  Alfven  waves. 

The  diagram  belov.-  indicates  the  various  subdivisions  of 
phenomena  vrhich  are  induced  according  to  the  above  scheme.   The 
superscript  in  parenthesis  outside  of  a  circle  indicates  the 

order  of  the  included  root,  e.g.,  WE)^~^      represents  the 

2     P  P  2 
Vacuum  ElectrodynamAc  roots   (o)  -  k~c~)   =0.   A  bar  over  a 

classification  symbol  means  that  dissipative  effects  are 

included.   The  star  in  I'lHD*  indicates  the  inclusion  of  the 

displacement  current.   The  s^nnbol   LR^"''^   is  indicative  of 

the  N  '  order  Lar-m.cr  resonance,  viz.,   cd '^  Nfi  +  ...   (IT  = 

integer,  .Q.  -   Iarm.or  frequency). 

The  first  relevant  observation  is  that  the  PEM  class  of 
phenomena  exhibits  a  temperature  decay  (TD)  to  aT  =  Z,    vrith 
a  decay  tim-e:  {2-y')~~.      The  superscripts  (O)  and  (2)  on  the 
TD  s^nnbol  indicate  that  in  the  ?£!•!  domain  the  displacement 
--■zv:.     iT  =  0  is  a  zeroth  order  effect  in  k  vrhile  in  the 
M-Z  l:main  it  is  a  second  order  effect  in  >.   If  vre  look  ahead 
to  the  m.atrix  form.ulation  of  the  problemi  for  a  mioment ,  then  it 
is  quite  evident  that  to  lov,"est  order  the  master  matrix  is 
diagonal  -.■rith  a  2x2  block  multiplying  the  tv:o  temperature 
entries  in  the  system  vector.   It  folloivs  that  in  "r/.e  ?Z::  3ase, 
the  system  eigenvector  "which  corresponds  to  the  os:illi7:ry  Zi 
modes  includes  1  -  ~  =   Z.      The  c om.pl ementary  eigenvector 
includes  a  tem.perature  decay  alor-g  the  path  T  4-  T  =  C,   i.e., 
the  perturbation  temperatures  are  disolaced  equally  and  oppositely 

The  J^jnction  betv^een  MHD  and  MHD*  indicates  that  the  effects 
of  the  displacement  current  becom.e  nesrlisible  for  v/aves  v/hose 
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velocities  are  much  less  than  the  speed  of  light,   c.   The 
formalism  subsequently  presented  Indicates  that  the  lowest 
order  form  of  the  current  equation  is  the  classical 
E  +  U  XB  =0.   If  these  roots  were  extended  to  next  order  in 
k,   both  the  dissipation,   -v-^^J  and  the  Hall  current  term, 

n  JXB   (cf.  Eq.  B-II,  23)  would  appear,  in  addition  to 

20  21  22 
the  latter  form.   However,  most  investigators   '      disregard 

the  JXB  term,  when  examining  these  disslpative  effects.   Such 

results  (i.e.,  excluding  the  Hall  current)  would  follow  from. 

the  more  general  theory  in  the  limit  of  large  number  density. 

The  second  region  of  phenomiena  is  governed  by  the  disper- 
sion relation  go  =  cd^*^^  +  kco^''"''  +  ...,  o)^"'''  ^   0.   As  already 
noted,  the  param.eter   ck  which  enters  from  the  Maxwell's 
equations  part  of  the  system,  is  treated  as  an  independent 
parameter,  which  accounts  for  the  presence  of  k   in  the  lowest 
order  dispersion  relations. 

Pure  oscillations  result  in  the  limit   ck  — ^  0,   together 
with  k  ^  0.   This  limit  gives  rise  to  decaying  plasm^a  o&cll- 
lations   (PC^'^'^),  and  decaying  plasma  magnetic  oscillations 
(PMO^^'').   These  roots  are  identical  to  those  found  by  Grad 
who  dem-onstrates  directly  that  the  plasma  oscillations  are 
associated  with  current  fluctuations  parallel  to  B^,   while 
the  plasm.a  magnetic  oscillations  are  associated  with  current 
fluctuations  perpendicular  to  3^.   In  general,  in  the  following 
limits  wherein  B  ^0,   those  modes  which  do  not  contain  the 
magnetic  field  are  the  eigenvalues  of  the  eigenvector  (J^^,  O), 
while  the  complementary  non-degenerate  roots  are  the  eigen- 
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values  of   (O,  Jj_),  where   J   Is  the  current,  and  directions 

4 
are  with  respect  to  B  .   The  remaining  fourth  order  cd  =  0, 

Is  the  degeneracy  of  the  vacuum  electrodynamlc  roots 

{o^    -   k^c^)^  =  0,  in  the  above  limit. 

The  problem  of  Larmor  resonances  has  been  estenslvely 

6,7,24 
covered  in  the  literature.        However,  most  authors,  when 

Investigating  this  problem,  impose  the  constraint  of  stationary 

ions  and  do  not  include  the  effect  of  collisions.   In  this  vein. 


investigated.   The  lowest  order  non-trivial  dispersion  of  the 
N   order  resonance  includes  terms  of  order  \^^   ~    .   The  effect 
of  the  self-collisions  alters  the  zeroth  order  (k  =  O)  results, 
wherein  the  decay  time  is  as  the  Inverse  self-collision 
frequency  v,   and  is  independent  of  the  resonance  order  of  the 
mode.   More  generally,  the  dispersion  was  seen  to  vanish  in  the 
limit  of  large  magnetic  field  or  vanishing  number  density. 
These  Larmor  resonances  are  seen  to  constitute  a  doubly  degen- 
erate four  fold  infinity  of  roots,  two  about  every  multiple 
Larmor  frequency  of  both  species  (cf.  footnote  p.  50). 

In  the  more  general  case   kc  =  constant,   k  =  0,   one 
obtains  the  lowest  order  propagating  PEM  phenomena.   The 
governing  equations  for  this  region  include  a  full  equation  of 
motion  for  the  current  as  opposed  to  the  simplified  E  +  UXB  =  0 
or  MHD.   However,  as  already  observed,  in  the  limit  of  large 
number  density  (i.e.,  plasma  frequency  large),  the  m.ore  general 
equation  of  motion  for  the  current  collapses  to  the  above  form. 
This  again  gives  substance  to  the  idea  that  the  Lund qui st 
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equations  govern  a  fluid  dynamical  medium,  since  the  relevant 
equations  follow  from  the  more  general  ones  In  the  limit  of 
large  number  densities.   In  a  similar  manner,  the  more  general 
momentum  equation  of  MHD  (cf.  Eq.  C-IV,  12)  collapses  to  the 
simpler  momentum  equation  of  PEM  (cf.  Eq.  C-IV,  29)  in  the 
limit  of  vanishing  temperature.   In  these  connecting  limits 

the  sound  speed  is  obviously  surpressed,  so  that  the  fast 

2     '2  / 
wave  of  MHD  becomes  the  Isotropic  wave,   v   =  A    (cold  fast 

wave),  while  the  Alfven  waves  remain  intact.   The  remaining 

roots  in  this  limit  of  large  number  density  are  the  plasma 

oscillations  (PO^  ')    and  four  plasma  magnetic  oscillations 

(PMO  ),  (cf.  Eq.  C-IV,  56).   If  the  plasma  frequency  is 

increased  far  beyond  the  Larmor  frequency,  then  these  roots 

2 
to  to  a  degenerate   o)  =  0,   together  with  the  two  plasma 

2    2 
oscillations,   cu  =  03  .   In  the  opposite  extreme  of  the  Larmor 

frequencies  far  exceeding  the  plasma  frequency,  these  plasma 

magnetic  oscillations  reduce  to  the  four  pure  Larmor  modes 

(lr(^'). 

The  effect  of  decreasing  the  large  number  density  in  this 
region  was  uncovered  by  extending  the  related  calculation  to 
terms  of  order  [kc/co  ]  .   In  the  direction  parallel  to  the 
magnetic  field  the  speed  of  the  cold  fast  wave  becomes  separ- 
ated from  that  of  the  Alfven  wave  by  the  amount 


(kc/co  )((m-m)//m  m) ,      where   m  and  m  are  the  ion  and  electron 
masses  respectively.   In  the  normal  direction  the  velocity  of 
the  cold  fast  wave  is  decreased  by  the  factor  (l-(k  c  /go  ) ) . 
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The  first  directed  wave  (kx.B  =■•  O)  yields  eight  decaying 
plasma  magnetic  waves  (PMW^^)  ^  ^    (cf.  Eq.  G-IV,  64)  and  two 
decaying  plasma  oscillations.   These  waves  are  the  extension  of 
the  classical  non-dlsslpative  waves  reported  by  Bernstein  and 
Trehan.    In  the  limit  of  large  number  density  the  magnetic 
field  both  the  cold  fast  waves  and  Alfven  waves  are  recaptured. 
Spltzer   has  discussed  these  roots  in  the  limit  of  large  mass 
ratio  in  which  case  it  is  impossible  to  recapture  the  correct 
MHD  roots.   This  large  mass  ratio  limit  is  also  reported  by 
Oster'  who  recalls  that  they  are  the  ordinary  (+)  and  extra- 
ordinary (-)  modes  of  ionosphere  theory. 

The  second  directed  wave  (k-B^  =  O)  (cf.  Eq.  C-IV,  69) 
again  gives  the  extension  of  the  non-dlsslpative  waves  reported 
by  Bernstein  and  Trehan' ' .   The  six  plasma -magnetic  waves 
(PMW_j_)  '      in  the  limit  of  large  number  density  and  large 
magnetic  field  yield  the  cold  fast  waves  in  the  direction 

9  =  7r/2.   The  decaying  magneto-ionic  roots  (MI^  '^)   are  a 

2S 
generalization  of  the  non-dlsslpative  magneto-ionic  roots.  -^ 

The  root  that  degenerates  to  co  =  0  when  the  dissipation  is 

2  2    2     2  2 
set  equal  to  zero  behaves  as:   Icd  ~  v-,  c  k  /  (co  +  c  k  )   (cf. 

Eq.  C-IV,  44.1),  for  small  collision  frequency   V-,  . 


Cf.  Ref.  20,  Eq.  III-7I 

r 


Cf.  Ref.  16,  Eq,  4-13 


Cf.  Ref.  25,  Eq.  131 . 
TT   Cr.  Ref.  20,  Eq.  III-82. 
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For  waves  propagating  In  an  arbitrary  direction,  two 
limiting  cases  were  considered.    These  are  vanishing  number 
density  and,  secondly,  vanishing  magnetic  field.   In  the  first 

limit  (cf.  Eq.  C-IV,  39), as  expected,  the  vacuum  electromagnetic 

(h)  2 

roots  (VE^  ')    are  recaptured.   Besides  a  degenerate   co  =  0, 

the  remaining  modes  are  a  purely  decaying  mode  and  four  decaying 

Larmor  modes  (LR^  ').         (The  superscript  (l)  Is  Indicative  of 

first  order  resonance.) 

In  the  remaining  limit  of  vanishing  magnetic  field  the 

decaying  plasma  oscillations  and  decaying  magneto-ionic  waves 

2 
once  more  appear,  in  addition  to  the  doubly  degenerate  cd  =  0. 

The  analysis  which  is  to  follow  suggests  a  connection 

between  the  so  called  macroscopic  and  microscopic  approaches. 

The  distinction  between  these  two  regions  is  that  the  former 

domain  is  governed  by  the  distribution  functions,  while  the 

latter  is  characterized  by  the  less  detailed  moments  of  the 

distributions.   In  the  present  formulation  the  distribution 

function  appears  as  a  linear  functional  of  the  first  three 

moments  of  both  gases.   This  functional  is  expanded  about 

small   k.   The  m.ore  terms  that  are  retained,  the  more  precise 

is  the  resulting  approximation  to  the  distributions.   At  any 

point  in  the  procedure  macroscopic  equations  may  be  obtained 

by  forming  the  proper  inversions  back  to  (x,t)  space.   The 

order   of  these  differential  equations  relates  to  the  order 


Since  there  are  only  three  external  parameters  in  the 
problem,   Bq,  Uq,    TqJ   and  Tq   is  already  constrained,  there 
remain  only  two  relevant  extrem.e  possibilities. 
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of  k  kept  In  the  expansions.   In  this  sense  the  detailed 
fluctuations  of  the  distributions  are  characterized  by  the 
higher  order  derivatives  in  the  macroscopic  equations.  A  good 
example  of  this  connection  is  in  the  magnetic  resonance  problem, 
wherein  the  dispersion  of  the  N'^^  resonance  is  of  order  k 
so  that  the  macroscopic  equations  of  order  2N-2   exhibit 
resonances  up  to  order  N. 


B.   Equations  and  Linearization. 

I .   Linearization  of  variables. 

Before  wading  into  the  formulation,  perhaps  it  is  wise 
at  this  point  to  impart  a  note  on  the  notation.   To  distinguish 
between  ion  quantities  and  electron  quantities  we  adopt  the 
convention  that  ion  variables  carry  a  roof  and  electron 
variables  are  bare.   For  instance,  the  two  macroscopic  veloc- 
ities are  written  u  and  u;   the  former  for  ions  and  the 
letter  for  electrons.   To  distinguish  between  full  variables 
and  perturbation  quantities  we  put  a  bar  on  top  of  the  full 
variable  and  nothing  on  the  perturbed  variable.   For  instance, 
the  perturbation  of  the  electron  temperature  is  described  by 
the  equation  T  =  T  +  T.   All  variables  that  make  up  the 
initial  equilibrium  situation  carry  a  subscript  zero  which 
explains  T   in  the  latter  equation. 
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The  linearization  is  about  an  absolute  Maxwellian   f^  ' 
and  obeys  the  following  scheme: 


(1) 


.(o) 


E  =  0+E,      n  =  n+n,   u  =  0+u,     p..=  P^*^  + 


ij    ^o   '  ^-^ij 
B  =  B   +B,     T  =  T+T,   U  =  0  +  U,     "C^Cq+C 
(2)   w  =  w   +  w,     p  =  p  +p, 

J  =  0  +  J,      Q=0+Q,    p=     p^+p,         N  =  N^+N 

P.  .  =  P-^'^  +  P.  .,         Z  =  Z+Z,    W  =  W  +W. 
ij    o     ij'         o    '         o 

Except  for  the  individual  temperatures   T,  T,   capital  letters 
denote  combined  fluid  variables.   They  are  related  to  the 
lower  case  letters  (individual  flow  fariables)  in  a  purely 
additive  manner,  i.e.,   A  =  a  +  a    (A   represents  p,    P,  N, 
Z,    W)   except  for  the  mean  flow  U  which  is  given  by: 


(5)  pU  =  pu  +  pu, 


In  the  above  set  (2),   E   is  electric  field  intensity, 
B   is  magnetic  field  intensity,   w   is  internal  energy,   n 
is  number  density,  p  is  mass  density,   u   is  macroscopic 
velocity,   p   is  pressure,  C,      is  heat  flow,   J   is  current 
density,  and   Q  is  charge  density.   It  should  be  noted  that 
to  any  individual  flow  variable,  (lower  case  letter)  which 
does  not  carry  a  roof  (i.e.,  electron  variables)  there  corre- 
sponds a  similar  symbol  which  is  roofed  (ion  variable).   An 
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example  is  u  which  is  the  said  "reflected"  u  variable.   In 
this  manner  there  are  eight  more  equations  in  the  set  (1,2) 
which  are  so  generated. 

The  equilibrium  configuration  is  further  defined  by: 


N  M 
^0+  "^o  =  ^ 


w  =-JnkT,    w  -Jn  kT^,   T  =  T^,     so  that 
o   2   o   o'     o   2   o   o    o    o 


1 
'o  '  '"o    2  ^'o'^'o    "o 


w.  +  iAf_  =  ^  N_kT,  =  W_  ,     w   =^  w^  =  wy: 


It  is  convenient  at  this  point  to  introduce  the  scalar 
pressure   p. 


(5)  p=iTrp..  =|w. 


The  notation  is  such  that  any  p  without  subscripts  is  a 
scalar  pressure. 

This  latter  identity  symbol  follows  from  the  kinetic 
definition  of  pressure  and  internal  energy  (to  be  introduced 
momentarily) . 

The  stress   s.  .   is  defined  as  the  traceless  part  of  the 
pressure  tensor,  viz., 

(6)  s.  .  =.  F^j  -  P  . 


For  the  equilibrium  scalar  pressure  we  write 
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Po  +  Po  =  ^o  '         Po  =  Po  =  V2  -  ^o^^o 
(7) 


The  absolute  Maxwellian  f^°'      is  about  zero  mean 
velocity  and  is  given  by 

(8)        A°)    ,-^^e.,{-e/2C^]^"-j^    , 

which  also  defines  the  exponsntial  function  j6.  The  particle 
velocity  is  i,  and  the  thermal  speed  C  is  defined  through 
the  equilibrium  temperature   T  ,   viz.. 


(9)       mC^  =  kT^  =  mC^. 


The  electron  and  ion  masses  are   m  and  m  respectively. 

Let  us  now  formulate  expressions  for  these  perturbation 
quantities  in  terms  of  velocity  moments  of  the  perturbation 
distribution  function  g.   Firstly,  for  the  number  density  we 
recall  the  kinetic  definition: 


(10)       n  =  /  f  d^l  =  /  (f ^°^  +  g)d^^  =  n^  +  n 


so  that 


(11)      n  =  /  g  dh. 


The  mass  density  is  related  to  the  number  density  as 
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N_M 


(12) 

(12' )     M  =  m  +  m, 

so  that  the  perturbation  mass  density  Is  given  by; 

(13)      P  =  nm  +  nm  . 

For  the  velocity  u  we  recall 


(14)      Hu  =  /  ifdh   =  /  ^(f^°^  +  g)d^4 


(15)  nu  =  0  +  n^u  , 
so  that, 

(16)  n^u  =  /  4gd^|  ,    and  p^u  =  m  /  ^gd^^ . 

It  follows  that  for  the  mass  averaged  flow  U,   we  have 

(17)  U  =  0  +  U 
and 

(18)  p^U  =  p^u  +  gG  . 

The  Internal  energy  Is  defined  In  terms  of  the  kinetic 
energy  measured  in  a  frame  which  is  moving  with  the  velocity 
U,   so  that 
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(19)      W  =  ^  m^.  /  fj4-U)V4 
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The  sum  over  repeated  indices  covers  Ions  and  electrons.   The 
linearization  yields: 

W  =  W^  +  W  =^  I  m.  /  f|°hV^  +  I  m.  /  g.^V^  , 


(20) 


(21)      w  =^  I  m  /  g^^dh    ,  w^  =  I  n^KT^ 


We  recall  that  the  temperature  Is  related  to  the  Internal 
energy  through  the  kinetic  definition  of  temperature  (internal 
energy  per  particle): 


(22)      w  -  ^  nKT 


the   linearization   of  which  gives 


(23) 


nKT   =  K(n     +   n)(T     +   T)    =   Kn   T     +    K(nT     +   n^T) 


im/e^f  =  i,n/e2(f(°)+g), 


whence  It  follows  that 

(24)      ^  m  /  l^g  =  k/  gT^  +  n^TK  . 

In  terms  of  the  thermal  speed   C   =  kT/m  this  latter  equation 
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appears  as 

n  TC^ 
(25)             -^  - 

0 

_   1 
3 

/  g(^^-  3c2)d^^  , 


the  non-dimensional  form  of  which  appears  as 


(26)      T  =  ^  J   g(|^-  3)d^^  . 


If  the  subscript  d  denotes  "dimensional",  then  T,  =  '^n'^' 

n 
g  =  — :t  g  ,   ^H  ~  ^/C  ■   We  will  return  to  the  non-dlmenslonal 

formulation  when  discussing  the  Boltzmann  equations. 

The  pressure  Is  defined  In  terms  of  momentum  transport 

measured  In  a  frame  moving  with  the  velocity  U.   We  write 

accordingly 


(27)      P,,  =  m.  /  f.(|-u)^(|-U)^ 
Linearization  yields 


m,  /  4°)  ^/,  +  m^  /  g,^,^. 


(28) 


Note 


(29) 


P   +  P„ 
o    In 


P^  =  i  Tr  m.  /  f|°)  e,^^  =  ^  Tr  P 


P^n  "  ^n  /  S,^,^^  =  P, 
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similarly,  the  heat  flow  vector  Z   Is  the  heat  flow  measured 
In  the  U   frame . 

(30)       Z  =.  |m.  /  f.(4-U)2(e-U)  . 

Linearization  yields  . 

Z  =  -U  I  m.  /  f|°)(^2  ^  j^2^    +  I  m.  /  g^iH 


(31) 


Z  =  0  -  U  I  (N^kT^)  +  I  m^  /  gi^^^ 


and  for  the  Individual  heat  flow  we  write: 
(52)  C  =  0  -  U  I  n^KT^  +  I  m  /  gl^^  . 
It  follows  that  the  perturbation  quantities  are  given  by 

(33)  z  =  -|uw^  +  |m.  /  g.|2^   , 

and   for  the   individual   perturbation  heat   flows: 

(34)  c  -  -  |uw^  +  |/  mgi^i     . 

It  is  proper  at  this  point  to  introduce  a  less  physically 
significant  heat  flow   C ' ^  the  energy  transport  measured  in 
the  individual  frames  of  each  gas.   This  heat  flow  is  generated 
from  the  above  formulation  by  letting  U  — >  u   for  each  compon- 
ent, so  that   C'   appears  as 
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(35) 


^  '  =  -  I  uw^  +  I  m  /  gl^l 


C'  =  |/  g(^^-  50)Uh. 


It  follows  that  the  difference  between   C   and   ^'   is 

w  m 

m-f-m 
(36) 

C  =  C'  +  ::^  (u-u) 

m+m 
so  that  the  difference  In  the  total  heat  flows  Is 

(37)  Z  =  Z'  -  w^   1^   (G-u)  . 

m+m 

These  two  heat  flows  will  be  the  same  if  either  the  masses  of 
the  two  components  are  equal  or  the  velocities  of  the  components 
are  equal . 

In  the  non-dimensional  units  introduced  above,   C'   appears 
as 

(38)  C  =  I  mn^C^  /  g(|2_  ^)^^\    . 

The  charge  density  Q  is  given  by: 

Q  =  e(n-n) 
(39) 

Q  =  0  +  e(n-n), 

so  that  the  perturbation  charge  density  is  written 

(^0)       Q  =  e(n-n). 
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The  current  density  J   is  written 
J  =  e(nu  -  nu) 


(^1) 


J  =  0  +  en  (u-u) 
o  ^    ' 


whence  the  perturbation  is  given  by 


(42)      J  =:  en  (u-u) 


Before  leaving  this  section  on  the  formulation  of  the 
perturbation  variables,  we  note  for  future  reference,  the 
following  relation  between  the  individual  fluid  variables 
[n,n;  u,u],   and  the  total  fluid  variables   [  j-^,Q;  U,J]. 
First  we  repeat  the  above  derived  relations. 

j^=mn  +  mn  J  =  en  (u-u) 

(^3) 

Q  -  e(n-n) 


and  secondly,  the  inverses 
(W) 


u    = 

=  -e„:«^ 

^ 

u    = 

TT                   m         -r 
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B-II.   Boltzmann  and  Moment  Equations. 

In  Appendix  A  we  Indicate  the  manner  In  which  the  pertur- 
bation away  from  the  above  said  equilibrium  state,  i.e., 
[f,f,E,B]/   ^^-^N  =  [f^°S  f^°  ,  0,  B^],   and  subsequent 
linearization  yields  the  following  linearized  Boltzmann  equa- 
tions . 

(1)  vg  +  |£.e.vg-^exB„.v5g-|E.V5f(°) 

n=0  ™  ^ 


oz  m      "   ^    m     ^ 

n=:0  m 


The  H^  ^   Hermlte  polynomials  are  generated  from  a  local 

13  ih 

Maxwelllan   '    In  the  following  manner.   If  the  local 
Maxwell Ian  f 


(3)       f  - ^^7o  exp  [ 


(27rRTp^         2RT 


Is  expanded  about  the  equilibrium  configuration,  linearization 
yields 

(^)        f,  =  f(°)[l  +^h(°)  +  ^  •  H^l)  +  ^h(2)] 
°  ^o        ^  ^o 

where 

(5)       H°=l;   H(l)=i;    h'^)  =  |  [ (|)2-  3]  . 
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H^^'^  accordingly  Is  a  function  of  £,/C.   The  coefficients  a^'^'^ 
are:   a^^''  =  n/n^;   a^^^  =  u/C;   a^^^  =  T/T^. 

The  variable ''<-  is  the  reduced  mass  mrn/(iii+m).   We  note 
that  the  model  includes  four  dissipation  constants,  viz., 
v,v,v-,,v  .   The  first  two  are  related  to  the  self -collision 
frequency  of  the  Krook-Bhatnager-Gross   model  and  are 
associated  with  the  time  of  relaxation  of  the  gas  to  a  local 
Maxwellian  state.   The  frequency  v,  emerges  as  being  related 
to  the  electrical  conductivity,  and  is  associated  with  a 
relaxation  from  a  state  of  flow  of  one  gas  through  the  other 
to  a  quiescent  state.   The  third  dissipation  frequency   v^ 
Induces  an  equilibrium  state  wherein  both  gases  have  the  same 
temperature,  i.e.,  the  same  energy  per  particle. 

These  constants  are  positive  so  as  to  Insure  relaxation 
in  the  correct  direction  in  time.  In  addition,  the  specific 
orientation  of  the  mass  parameters  will  be  seen  to  yield  the 
correct  conservation  equations. 

Let  us  here  recall  the  integral  forms  of  the  flow  variables 
derived  above: 

n  =  /  gd^^  ,    n^u  =  /  egd^l  ,    w  =  I  /  g^2^3^ 
(6)        n^(T/T^)  =  \  !   g((i)'-  ^)ah    ,    P,,  =  m  /  g^.^^d^^ 


C  =  §  /  g(^^-  5C^)idh    , 


which  also  serves  to  generate  orthogonal  polynomials  with 
respect  to  the  Maxwellian  weighting  factor  f^  '^.   To  preserve 
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orthogonality  of  these  polynomials,  the  temperature  is  taken 
as  the  third  moment  Instead  of  Internal  energy. 

Returning  to  equations  (l),  we  consider  the  three  opera- 
tions: J  (P^,      J   d^U^,   /  d^i  f  i^ ,      to  obtain 

(a)   The  continuity  equations: 


(7)       If  +  n„  Vu  =  0 


(8)       ||H.n„VG  =  0. 


Where  we  have  recalled  that  V^f^°'^  =  -^f^°VC^   and 


(9)       IXB-  V^g  =\7^.exBg 


(b)   The  momentum  equations: 


When  we  operate  with   /  £, .  d  |   there  remains 

e  "^1 

(11)       1^  /  dh    l.^.g  =  1^  ^     by  Eq.  (B-I:  29). 


The  change  Is  one  of  normalization.   Instead  of  Including 
the  energy  per  unit  volum.e,  w,  we  discuss  T,  the  energy 
per  unit  volum.e,  per  particle. 
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The  ^^B     term  Is  expanded  as  follows: 

which  Is  the  1  ^  component  of   2,  X;  B,   so  that  the  above  moment 
equation  reduces  to  (here  we  Include  the  ion  equation) 

(12)      n^'^^l   njE  +  .xBj  H-  i  Vp„  = 'f  v^n^Au 

Let  us  observe  in  passing  that  the  formation  of  these 
moment  equations  supposes  the  vanishing  of  g  at  the 
boundaries  of  velocity  space.   It  will  be  interesting  to 
observe  the  location  of  this  assumption  in  the  analysis  section 
to  follow.   One  also  notices  that  the   Au   interaction  terms 
indicate  that  the  relative  motion  of  the  gases  serves  as  a 
momentum  source  for  the  individual  components. 

(c )   The  energy  equations: 

To  form  the  transfer  equation  for  w  we  operate  with 

1       ^2 

^  m  /  d  2,?,  .   First  let  us  note  the  manner  in  which  the  self- 
collision  terms  cancel. 

(14)     i  vm  /  f^°k^{^   +  ^  -<-  ^  H^^bA  =  V  I  k(T^n  +  Tn^) 
o   C      o     . 

=  vw  =:  I  mv  /  ge^d^^   , 
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so  that  the   v   terms  cancel,  and  the  only  remaining  collision 
term  Is  the  AT  term. 

The  magnetic  term  Includes  the  factor 


(15)        /e'v-exBs  =  /?2|^e„B„B^„,g 


=  !  kj   (^^^nVnm^s)  "  ^SmA  ^   S^n^r 


The  first  term  vanishes  because  It  Is  a  surface  term,  while 
for  the  second  term  we  note  that  the  vector  representation 
Includes  the  factor  ^X  B-^,      which  vanishes  since  ^XB 
Is  normal  to   I. 

Finally  for  the  spaclal  part  of  the  convectlve  derivative 
we  have,  by  equation  (B-I:  3^) 


(16)      V.  |m  /  g^h   -   V-iZ  +   |uw^)  . 


The  energy  equations  appear  as 


(17)      It  ^  +  f  ^0^°"  +^°^  ^  i  ^2   ^o^^^ 


(18)      It  ^  +  f  w^  V-U  +  V"?  -  -  I  v^  n^kAT  . 


Here  we  notice  that  a  displacement  In  temperature  between 
the  two  gases  acts  as  an  energy  source  for  the  Individual  com- 
ponents. 

It  Is  consistent  at  this  point  to  form  some  linear  combina- 
tions of  the  above  equations  (7,8;  12,13;  17jl8)  so  as  to  obtain 
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a  more  familiar  set.   For  the  most  part,  we  are  making  the 
transformation  of  variables  (cf.  Eq.  B-I:  43,^4). 

n,  n  -^  P,  Q 
(19) 

u,  u  -^  U,  J  . 

Mathematically,  either  system  has  its  independent  desirability, 
but  physically  the  transformation  is  significant  because  any 
observations  which  are  made  in  the  laboratory  are,  indeed, 
with  respect  to  the  variables  on  the  right.   We  note  that  the 
list  is  not  extended  to  the  higher  moments  (e.g.,  p..,  C>  .-•)• 
Although  one  readily  agrees  that  the  total  components  are 
physically  relevant,  e.g.,  P.  .  =  p.  .  +  p..  .,   it  would  appear 
that  the  related  difference  variables  (with  appropriate  mass 
factors)  have  not  come  to  the  fore,  so  to  speak.   These  combina- 
tions will  suggest  themselves  immediately  when  we  form  the  so 
called  difference  equations. 

To  obtain  the  system  equivalent  to  (7^8)  we  multiply  first 
by  the  respective  masses  and  add.   The  second  equation  is 
formed  by  multiplying  by  the  charge   e   and  subtracting  the  ion 
equation  from  the  electron  equation.   In  this  manner  we  obtain: 


(20)       l-p  +  P  V-U  =  0 


(21)       1^  Q  +\7-J  =  0. 


In  the  exact  same  manner  we  form  the  system  equivalent  to 
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the   two  momentum  equations: 


(22)  poltU+    V-P 


ij 


(23)  U  +   \7  -R.  .    -  00^    £    (E  +  U  X  B    )    +  n   J  X  B  =    -V,  J 

^    "^ '  dti.ioo^  o'  1 


To  obtain  the  second  equation  we  noted  that 


(24-)      u/m  +  u/m  =  U/^  +  -^   [m-m/m  m] 

o 


(25)      ^  =  m  V(m+m)  . 

In  the  second  of  these  equations 

2 
constant,   cd   Is  the  plasma  frequency, 

measured  with  respect  to  the  reduced  mass  9^1  .        fi   Is  a 
modified  Larmor  frequency. 


(26)      fi  =  eB^(m-m)/m  m  =  fi-fi 


and  B  Is  a  unit  vector  In  the  B  direction.   Lastly,  the 
difference  quantity  R.  .   Is  given  by 


In  a  more  Inclusive  formulation  (really  notation)  the  term 
(+e)  should  be  replaced  by  e,  to  account  for  multiply 
charged  ions.   Uncer  such  circumstances  R. .   would  keep 
the  above  form  with  e/m  replaced  by  e/m.   "''"' 
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We  note  that  In  transforming  from  the  two  Individual  momentum 
equations  to  the  above  two  equations,  (so  called.  Momentum  and 
Ohm's  law  (or  current)  equations  respectively),  the  Individual 
pressures  were  transformed  according  to 

(28)        P.J,  5..  ^p^.,  R.^.  . 

This  appears  to  be  a  natural  extension  of  the  transformation 
set  (19)  In  so  far  as  one  of  the  two  new  variables,  viz.,   R.  ., 
is  of  an  inherent  electric  nature.   A  direct  interpretation  of 
R  is  obtained  in  the  limit  where  Ohm's  law  reduces  the 
conservation-like  form 


(29)       |^J+7-R,j=0, 


SO  that  the  time  Increase  of  the  current  density  is  associated 
with  a  negative  divergence  of  R.  .,   or  a  net  influx  of  the 
vector   R.  .  ,  (into  the  related  volume  element).   This  is  not 
surprising  if  we  recall  that  the  pressure  tensor  p .  .   is  the 
rate  of  transport  of  i-momentum  (of  the  ions)  in  the  j-direction. 
If  we  multiply  by  e/m  we  obtain  the  rate  of  transport  of 
1-current  (that  part  due  to  the  ions)  in  the  j-dlrection. 
Subtracting  (e/m)p.  .  gives  the  tensor  R.  .   which  is  seen  to 
be  the  rate  of  transport  of  the  total  i-current  in  the  j- 
dlrectlon.   We  accordihgly  refer  to  the   R.  .   tensor  as  the 
J-flux  tensor. 
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Finally,  we  consider  the  addition  of  the  two  energy 
equations  (17,  l8)  and  obtain 


(30)     |_  w  +  I  w^  V-u  +  V-z  =  0. 


If  the  total  heat  flow  Z  =  C  +  C   vanishes,  then  this 
equation,  together  with  the  continuity  equation  (for  P  ), 
yields  the  linearized  form  of  the  adlabatlc  law: 


(31)     P/Po  =  I  P/p, 


We  have  recalled  that  ^P  =^W  =  w  +  w,   and  that 

W  =  2w  . 
o     o 

The  physical  significance  of  the  electric  difference  terms 
is  not  well  defined  for  these  higher  moments  and  we  accordingly 
do  not  form  the  counterpart  to  the  above  energy  equation  (30). 
Nevertheless,  it  is  interesting  to  note  that  the  list  of  new 
variables,  viz.,  (P,Q);  (U,j);  (P,R).  .;   suggest  the  proper 
form  for  the  relevant  difference  equation.   Noting  that  the 
scheme  is  such  that  the  div.  of  a  successive  term  Is  coupled 
with  the  time  rate  of  change  of  a  previous  variable,  e.g., 
P+J^V-U;   PoU+  \7-P;   Q+  V-J;      J+V-R;   and  continuing 
to  the  energy  equation:   l/2  (Tr  P)  +^-Z,   it  is  quite  natural 
to  extend  the  list  to  Include  the  form,  (l/2)Tr  R  +\/'K,      where 
K  is  a  first  order  tensor  and  might  appropriately  be  labelled 


mass,   M/2,  viz.,   a^  =  [^P/Splg  =  |  [Pq/Po^  =  f  (^^c/f) 
It  is  formed  in  Appendix  F. 
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the  R-flux  vector. 

The  scheme  is  physically  relevant  in  as  much  as  each 
succeeding  equation  serves  to  determine  a  variable  which  is 
involved  in  a  previous  equation.   In  this  sense  once  the 
original  set  is  fixed  upon,  e.g.,  (Q,^)?  the  proper  remaining 
higher  moment  variables  are  defined. 

The  reason  why  we  have  taken  time  to  formulate  certain 
well-known  macroscopic  equations  is  three-fold.   Firstly,  we 
have  demonstrated  that  our  linearized  Boltzmann  equations  do, 
indeed,  give  rise  to  the  standard  linearized  form  of  the  con- 
servation equations,  (i.e.,  Eqs.  20-23).   Secondly,  it  is 
revealing  to  observe  the  manner  in  which  the  stated  collision 
forms  appear  in  the  flow  equations.   Lastly,  it  is  Important 
to  have  these  equations  on  hand  for  identification  purposes  in 
the  analysis  section  to  follow,  wherein  the  formulation  is  con- 
structed about  the  distribution  function  written  as  functional 
over  the  collision  terms.   The  emergent  equations  are  linear 
combinations  of  equations  which  are  similar  to  the  above 
forms.   Comparison  of  these  forms  enables  one  to  identify  the 
transport  coefficients  although,  in  this  analysis,  we  are  not 
specifically  interested  in  this  region  of  investigation,  but 
rather  in  the  characterization  of  propagation  disturbances. 

The  standard  sets  of  the  first  three  conservation  equations ^ 
together  with  Maxwell's  equations  which  give   E   in  terms  of 


H.  Grad  has  suggested  changing  variables  in  the^  Boltzmann 
equation  from  f ,f  to  P,G:   F  =  m?  +  mf ,   g  =  §f  -  ef 
(Private  communication.)   (Of.  Appendix  F). 
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U  and  J,   constitute  thirteen  equations  In  the  25  unknowns: 
(n,  u,  w,  p^.,    C;  n,  u,  w,  p^^,  t,    E) .   Note,  w  -  (l/2)Tr  p^y 
The  standard  macroscopic  technique  which  Is  used  to  close 
the  system  Involves  postulatlon  of  the  transport  coefficients. 
I.e.,  phenomenologlcal  statements  relating  higher  moments  to 
lower  moments.   So  It  Is  that  one  defines  (for  simple  gases 
say)  the  viscosity  r]:      ? .  .   =   t\   SU./Sx.  ;  and  the  heat  conduc- 
tion A:   ^.  =  A  St/Sx. .   Such  statements  for  the  above 
system  would  reduce  the  number  of  unknowns  to  thirteen  and  the 
system  would  be  determinant.   Such  phenomenologlcal  connection 
formulas  are  the  essence  of  the  so  called  macroscopic  approach. 

A  more  self-contained  technique  Is  that  of  the  microscopic 
approach  wherein  one  relies  only  upon  the  distribution  function 
In  order  to  describe  the  relevant  macroscopic  fluctuations  In 
the  fluid.   Indeed,  there  are  regions  of  phenomena  where 
certain  giant  fluctuations  In  the  distribution  function  Itself 
are  observable,  (e.g.,  Larmor  resonance,  cooperature  phenomenon 
In  general).   The  distinction  between  macroscopic  fluctuations 
and  those  just  mentioned  Is  that  In  the  latter  the  resonance 
regions  are  seen  to  be  regions  about  singularities  In  the 
distribution  function. 

In  the  following  section  we  shall  endeavor  to  employ  a 
microscopic  technique  In  order  to  uncover  the  possible  modes 
of  development  which  are  excited  by  the  Initial  perturbation. 
The  starting  point  will  be  the  said  linearized  Boltzmann 
equations . 
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C.   Analysis. 

I .   Formulation  of  Homogeneous  Equations. 

We  rewrite  the  equations  (B-II:  1,2)  In  the  more  compact 
notation 


(1)  vg  +  ||+  ^  Vg  -  ^  ^xB^-v^g  -  f^°^[ho  +  e-h^] 


where 


Tj;f  =    vT    -    v^AT 


o 


2 
u      ^   ,,      Au  ^   eE 


(3)  '^l-^-^^f  -i|.        H2=|[!|-51 


^  C        m  C         mC 


In  general,  to  obtain  an  Ion  form  from  an  electrom  form 
merely  roof  all  variables  and  let  -e  — > +e .   We  will  accord- 
ingly work  only  with  the  electron  equation  realizing  that  at 
any  given  point  in  the  analysis  the  corresponding  Ion  equation 
could  be  obtained  In  this  simple  manner.   The  above  system  Is 
seen  to  be  Invariant  under  the  transformation  F  — >  +   F; 
-e  -^  +e      where   F   Is  any  variable  that  Is  labeled  (sign  for 
V-,  ,  Vp  terms)  . 

* ' r2T 

Note  change  In  notation  H^  ^  — >  Hp  • 
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In  order  to  solve  for  the  distribution  function  in  terms 
of  a  functional  over  Its  moments  we  will  work  In  the  transform 
k,a3  space  rather  than  In  the  actual  x,t  space.   This  Is 
accomplished  through  the  introduction  of  the  Laplace -Fourier 
transform 

00   00 

^  '  O  -00 

Where  X  represents  all  dependent  variables  which  are 
Included  In  the  formulation.   The  transform  of  the  electron 
equation  appears  as: 

(5)       [v-lo.  +  Ik.^  +  Q   |^]g"  =  r^°hhl   +  ^h*). 


Let   B   be  a  unit  vector  In  the   B    direction.   It 

o 

follows  that 


(6)       |j_^Bx(|XB)   and   kj_  =  Bx(kXB) 


The  angle  a   Is  such  that 


Throughout,  we  tacitly  assume  the  presence  of  the  Initial  con- 
dition terms.   The  distinction  between  the  resulting  dispersion 
relations  Is  one  of  words  only.   In  the  Inhomogeneous  case 
encountered  when  the  Initial  conditions  are  explicit,  the 
related  roots  are  associated  with  singularities  of  the  trans- 
form inversion.   Those  encountered  In  the  homogeneous  case 
insure  solubility  of  the  system.   The  roots  are  Identical; 
they  are  both  zeros  of  the  master  determinant. 

Very  simply  |XB^-\7^  is  the  projection  of  ^7 e  onto  ^XB^. 
In  cylinder  coordinates  with  i  parallel  to  B  ,  ^/'da  is 
the  component  of  Vp      which  is  normal  to  t      and  B^. 
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It  is  the  angular  cylindrical  coordinate  of  the  vector  ?, . 
The  electron  Larmor  frequency  n  is  given  by: 

(8)       n  =  -eBym  . 

The  related  ion  frequency  is,   fi  =  +  eB  /m  . 

Dropping  the  star  notation,  the  equation  is  further 
condensed  to  the  form 


(9)       [a  +  ik^^j_  cos  a  +  |^]g  -  ^  h(a: 


where   g   is  now  the  nondimensional  distribution  measured  in 
units  of  ^^„/C   and  j6      is  the  exponential 

(10)      j6  =   ^7^  exp(-  ^) 


and 
(11) 


v-ico  +  ik  2, 


(12)      h  =  h^  +  h^-e   ; 

h   is  a  function  of  a   in  as  much  as  it  is  a  function  of   ^ 
Both  h  and  g  are  periodic  in  a,  so  that 

h(a)  =  h(a  +  2-w) 

(15) 

g(a)  =  g(a  +  2tt)    . 
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The  Integral  form  of  equation  (9)  appears  as: 

471 


(14)   K  =.  W^) 
^^^>        ^  4if~ 

e   -1 


exp<aa'   — ^ —  [ sin  a  -  sln(a+a ' ) ] |h(a+a  )da 


The  closed  system  of  algebraic  equations  Is  now  obtained  by 
forming  the  moments,  (g  and  i   nondlmenslonal ) , 


(15)       n  =  /  gd^l,    u  =  /  g^d^^,    T  =  I  /  gHgd^^ 


with  similar  forms  for  the  Ion  terms. 

It  Is  convenient  to  put  the  above  Integral  Into  a  more 
compact  form.  This  Is  accomplished  through  Introduction  of 
the  rotation  operator  R^,   which  Is  defined  by 

(16)  R0f(a)  =  f(a  +  e)    . 

A  relevant  example  follows  from  the  Identity 

(17)  ^  -  B(B-e)  +  BX(^  XB)  , 
so  that 

Rgl  =  B(B-^)  +  sin  e(^XB)  +  cos  9   Bx(^XB) 
(18) 

R  ^1  =  B(B-e)  -  sin  0(^XB)  +  cos  0  BX(sXB). 


This  Is  most  simply  obtained  via  the  use  of  the  Integrating 
factor:   exp[aa  +  lkj_^j_  sin  a]. 
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Returning  to  the  algebraic  formulation  we  note  that  any  of 
the  above  three  moments  are  of  the  form 

(19)      <H.>  =  /  gH.(^)d^^  =  /  d^^   4J   J  da'expjaa' 

o 


Transforming  variables,   4  — >  4'  =  R   i/p^'   so  that 
?,  =  ^_ctr/2^''   w^  obtain,  (dropping  the  prime  on  a) 

^-       f 

da  exp^  aa  -  ^=p  [R  ^  a  ^ 


o         ^  2   2 


(20)         <H.>  =  /  d^e-  Tifrj  ^"  ^^p| 

o         ^ 

2   2     J 

If  we  further  set  it  =   a/2,   there  results  (dropping  the 
prime  on   £, ) 


3* 


(21)       <H,>  -  2  I  d''^  \^^y       I   d-^  exp^  25^ 


e  ^^^-1  ^ 
o 


27r 


21k^-^j_ 

^^ sm  ^[  h(R^^)H.(R_^^) 


At  this  point  an  Interesting  effect  of  the  presence  of 
the  steady  magnetic  field  becomes  evident.   If  we  recall  that 
iA   Is  a  cylinder  coordinate  of  the  vdctor  4,   then  the 
Integration  Involved  In  the  above  equation  Is  easily  interpreted 


^5 


as  being  over  a  four  dimensional  velocity  space.   In  this 
sense  the  magnetic  field  imposes  an  extra  weighted  velocity 
integration  to  the  three  dimensional  averaging  involved  in 
the  formation  of  moments  without  an  impressed  magnetic  field. 
Equations  (21)  with  H.  set  equal  to  1,  ^,  H^ ,  together  with 
Maxwell's  equations  constitute  a  closed  set  of  thirteen 
algebraic  equations  in  the  thirteen  unknowns (n,  u,  T;  n,  u, 
T,  E). 

Symbolically  this  system  may  be  written 


(22)       A.  .  X.  =  0  . 


As  indicated  in  the  introduction  we  have  chosen  to  call  the 
13  X  13  matrix  A. .,   the  master  matrix  of  this  problem  and 
the  related  determinant  the  master  determinant.   The  relevant 
modes  of  excitation  are  obtained  by  setting  det | A .  .  |  =  0. 
The  integral  coefficients  of  A. .   are  functions  of 

(03,  k,  ^,   n,  o)  ,  V,  V,    V-.  ,    Vp,  n  ,  T  ,  K,  m,  m,  e) ,   where 

2    2/^  * 

■f^  =   ck,   CO  =  e  n  /e  ??2.  ;   are  the  only  new  variables.   These 

latter  two  variables  enter  from  the  Maxwell  equations  part 

of  the  system,  which  are  written  as  (all  quantities  are 

perturbation  quantities,  recalling  that   Q,  =  J  =  O), 


CD   is  written  instead  of   e  .   Also,  the  constant  \i^ 

o  °     2 

appears  through  -p  =  ck,   since   c   ==  l/e  u  .  Q,      is 

written  for  B  . 
o 
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V-B  -  0  ,  V-e^E  -  Q 

(23) 

VXE  =  -B,  VkB  =  [x^J  +  E/c^ 


c   Is  the  speed  of  light  and   e   and   I-l   are  the  dielectric 
constant  and  permeability  of  free  space,  respectively. 

The  latter  three  equations  are  combined  to  give  the 
standard  wave  equation  for   e  E  , 


(2^)       c^V^e^E  -  £^E  =  j  +  c^^Q 


The  transform  of  this  equation  Is  (where  variables  are  now 
transforms) , 


(25)       .„E  =  i  ^if^ 
CO   -  -p 


If  this  equation  Is  combined  with  the  transform  of  the 
continuity  equation,  viz., 

(26)      ooQ  =  k-J  , 


one  obtains, 

^'^^      ^o^-looj   ^2_^2 

Recalling  that   J  =  en  Au,   this  latter  form  may  be 
rewritten  as 


"o  JcD^Au  --Pi^-Au)  i 
FTZd)     2     2      f 


(28)       E 
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If  this  expression  for  E   Is  substituted  Into  the 
inhomogeneous  vector  h-^  (cf.  Eq.  (3)),  there  results 


""o  /-/>(-/'-Au)  1 


(29)       Ch,  = 


The  frequency  v..  Is  given  by 

2 

/Wi  1<^   CD 


^^°)       ^  =  m   'l-':?7^2 


CD 


The  corresponding  Ion  term  Is  given  by: 


(31) 

7n 

103^    03 

o?-f" 

where 

%  =   e^^o/^o"^      • 

In  the  following  section  we  will  Introduce  f\     which  Is 
the  non-dlmenslonal  form  of  h,   and  Includes  frequencies 
measured  In  units  of  the  Larmor  frequency  fi;   similarly,  the 
frequencies  of  ti  are  non-dlmenslonallzed  through  fi.   This 
process  of  non-dlmenslonallzatlon  generates  a  convenient 
variable   V^  of  the  two-component  plasma  which  Is  a  function 
of  the  reduced  mass  z*^: 

(32)       -V^  =  (vVfi)  =.  {vJ^)    =     1  .    o 


3  -  ''^-'    -    '-y-'    -        n*^   0.2-^2 


Where, 

(33)      f^*  =  eB^Tn 
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(34)      4^^  =  1 0.2/^1  .  ^2/^  =  en/B^B^ 


The  frequency  Q.        is  the  Larmor  frequency  of  a  particle  with 
charge   e   and  mass  /?7  -   mni/m+m  . 

Expressing  h-.   In  this  manner  as  a  function  only  of 
u,u,   eliminates  E   from  the  problem  and  reduces  the  master 
determinant  to  a   10  x  10.   Indeed,  the  homogeneous  set 
A. .X.  =  0   Is  now  purely  represented  by  the  equation  (21), 
and  the  related  explicit  coefficients  of  A.  .   are  produced 
by  performing  the  Indicated  Integration  over  the  four- 
dimensional  velocity  space,  with  volume  element  d   i   di/ . 


C-II.   Magnetic  Singularities. 
Let  us  recall  equation  (C-I:  21) 


(1)        I_g  =  <Hp  =  2 


d^^   l^^l    0  dii   exp  I  2af 
21k^-^^sin  ifj  h{R^i)E^{R_^i) 


The  variable  ?i   represents  the  non-dimensional   h,   and 
includes  frequencies  measured  in  units  of  the  Larmor  frequency 
fi,   and  velocities  measured  in  units  of  the  thermal  speed  C. 
Similarly,   h   includes  frequencies  measured  in  units  of  the 
ion  Larmor  frequency  fi,   and  the  ion  thermal  speed  C.   The 
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kj_"|j_   term  Is  non-dimensional  and  has  been  written  for 
[kj_'^j_/f2]  .   In  this  term  k_L   Is  measured  In  units  of  the 
Inverse  Larmor  radius  L"  .   Note,  W  =  C .      Again,  i^     is 
non-dimensional  through  the  thermal  speed  C.   In  general,  in 
the  ion  equation  all  explicit  speeds  are  non-dimenslonalized 
through  the  ion  thermal  speed  and  likewise  all  speeds  in  the 
electron  equation  are  non-dimenslonalized  through  the  electron 
thermal  speed.   A  similar  statement  holds  for  non- 
dlmensionallzlng  explicit  frequencies  through  the  Larmor 
frequencies  Q,     and  H,  respectively. 

It  is  quite  clear  that  in  the  special  case  k^^  =0  (a  -^  a) 
I.   Includes  regions  wherein  the  related  integrations  are  not 

■X- 

well  defined.   Since  these  resonances  occur  about   a  =  IN, 
i.e.,  -ico  =  INfi  -  v,   we  see  that  the  related  root  is  of 
the  form,   co  =  od^j       +   '^iq^^   +  .  .  .  .   We  will  show  in  this  section 
that   m  =  2N-2,   for  N   order  resonance,  and  derive  an  explicit 


The  variable  r\^     is  defined  as: 


(1-)      a  =  IN  +  Tij^ 


We  will  investigate  the  limit  of  f]^  -^   0,   kj_  — >  0. 
Knowledge  of  ■nvr(k)   determines  the  dispersion  relation  about 
the   N    order  electron  Larmor  frequency.   (By  N   "order" 


Actually,  they  occur  about   a  =  t  iN;   however,  we  will 
restrict  the  analysis  to  a  =  +  IN,   N  >  0,   realizing 
that  an  image  formulation  exists  for  a  =  -IN. 
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we  mean  the  N   multiple  of  fi.)   It  should  be  noted  that 
in  the  analysis  to  follow  no  assumptions  are  imposed  per- 
taining to  the  mass  ratio,  so  that  the  investigation  has 
a  complete  image  formulation  about  an  ion  resonance. 
In  the  above  limit  the  integral   I.   becomes: 

(2)       ^^  ==  i?  /5^^^^  j>   d^(^N^  +  2if)expk2m 

+   2ik_,-^j_  sin  -^j  h(R^|)H^(R_^l)  . 

Specifically,  we  are  interested  in  the  velocity  moment  u  =^  I-,. 
We  shall  see  that  the  reasonance  phenomena  is  characterized  by 
only  those  velocity  moments,  owing  to  a  natural  splitting  of 
the  system  (i.e.,  the  master  matrix  becomes  diagonal  with  a 
6x6  matrix  multiplying  the  six-component  vector  u,u).   The 
most  elementary  case  of  such  fluctuations  would  occur  in  a 
cold  plasma,  with  a  constant-in-time  number  density,  in  which 
all  of  the  particles  were  rotating  in  the  same  direction  with 
say  03  =  Q.   In  this  case,  T  =  T  =  0,  n  =  const.,  n  =  const., 
but   u   is  rotating  with  angular  frequency  Q.      Note  that 
there  is  no  translation,  so  that  in  the  usual  sense   u   (or 
more  appropriately  j)  is  not  apparently  fluctuating.   Of 
course,  such  variation  in  J   is  detected  by  the  electro- 
magnetic  radiation  emitted.    Such  a  phenomenon  is  properly 
categorized  as  a  giant  microscopic  fluctuation,  even  though 


E.  G.  Harris,  personal  communication 
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the  formalism  suggests  a  macroscopic  variation.   Similar 
arguments  apply  to  pure  plasma  oscillations.   In  this  sense 
the  two  phenomena  are  related  In  as  much  as  the  plasma 
oscillations  are  associated  with  linear  vibrations  In  J, 
while  magnetic  resonances  are  associated  with  circular 
vibrations  in  J. 

The  electron  macroscopic  velocity  appears  as: 

(3)      "   =  ^     ^^  9  ^^^^N^   ^   2if]exp<2±m  +   21y   sin   1/  I 

[h^+  h^(B(B-e)  +  (sin  ^)(4XB)  +  (cost)Bx(4  XB)  )  ]  [B(B- 4  ) 

-    (sin   ^)(4XB)   +    (cos   ^)Bx(^^B)]     . 

The  variable   y  =  lj_"k_i_.   and  B  =  B  /|B  |  .   To  evaluate 
the  21/^  contribution  we  note  that  the  lowest  order  non- 
vanlshlng  terms  appear  when  kj_  =  0.   The  related  Integrals 
are  evaluated  by  rewriting  the  sins  and  cos's  In  their  exponen- 
tial forms,  and  observing  that   /  xe^^dx  =  [  e^Va^]  [xa-1  ]  . 
In  our  case,  over  the  2tt     cycle,   only  the  [xe^^/a]    term 
remains.   The  pertlnant  Integrals  are  listed  below: 


However,  magnetic  resonance  is  clearly  not  a  cooperative 
phenomenon,  but  rather  an  ordered  response  to  an  external 
influence . 
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(p  e^^^"^  if   sin  if  dif  =  27r/[^N^-   U 

(|)  e^^^"^  ^  cos   -^  d-^  -    -l47rN/[4N^-   1] 

(4)  (fi  e^^^'''  if   sln^if   dif   -   l7r/[2N(N^-    1] 

e^^^"^  if   Qos^if   dif   =    -l7r(2N^-    1)/2N(N^-    l] 

(j)  e^^^'''  if   cos  ^   sin  ^  d^   -   7r/2(N^-   1]    . 

The  first  non-vanlshlng  contributions  from  the  r]^ 
terms  are  obtained,  again  by  writing  the  exponential  form 
for  the  sin  and  cos  terms  and  recalling  the  integral  form 
for  Jj,:   (p  expilN-^  +  IZ  sin  if  \  dif   =  2ir   Jj^(-Z).   The  pertinent 
terms  are  again  listed. 


(  1  e 


i2N*.lZ  sin*  ^,„,, ,_,,„.[ j^^_^.j^_^^^ 


,|.^^12N*+1Z  smt^^^,,,  ,  -[J2N.1  -  J^N-ll 

(5)       f  e"^*+^^  =1"  *  sin^V-  dt  =  .  Jg,-  f[J,,,2  -   J2N-2i 

X      12Ni/^+lZ    sin   if      ^2.     ,  ,  ^    t      j.   ^   r  t  +    t 

Q)  e  cos  if  dif  =  TT  JoM+   o    UoMio  "•"   •J- 


2N      2    ■-    2N+2  2N-2' 


2NV^+iZ    sin   if 


The  argument  of  all  of  these  Bessel  functions  is 
Z  =  2kj_'2,j_.   (Recall   k  measured  in  units  of   L   ,   and 
i      in  units  of  C.)   The  expansion  of  J^^  appears  as 
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(5)       ..Mi)l\,,^  (-r^w^f" 


N 


N!         m=l     m! (N+m) J 


It  follows  that  as  k  — >  0,  the  lowest  order  terms  appear  In 
the  quadratic  terms 

(sln^T/^ 
cos2^       \     d^ 
sin  i/^  cos  i/^ 


~  2(2N-2)!  ^^i"^x)      <f  +1  !^  + 


while  for  the  2if   contribution  we  repeat 

''  sin^'^      \  /+1 

(8)       (f  e^iN^  ^   cos^^        d^  =   if -(2n2-  1)  1  . 

\  sin  i/   cos  i/'/  \  -iN 


Having  thus  performed  the  dif   part  of  the  integration  we 

now  form  to  the  d  £,  part,  and  immediately  note  that  in  this 

limit  the   h  B(B«^)  term  vanishes  because  the  related  integrand 

is  odd;  while  the  h  sin  '^     and  h  cos  1/     terms  are  discarded 
o  o 

because  they  are  higher  order  in  k  than  are  the  quadratic 
terms  in  sin  1/ ,    cos  i/ .      In  addition,  we  note  the  splitting 
between  perpendicular  and  parallel  (to  B  )  components.   Con- 
sider Uj_.   The  contribution  from  terms  like  h  *b(B"|) 
vanishes  because  it  is  odd.   Similarly  for  all  of  the  other 
cross  terms  in  i    i^_'      At  this  point  it  is  convenient  to 


5^ 


Introduce  the  ||f||  symbol  defined  by 

(9)  llm    1  (()e21NV'+iy  ^^^  ^  F(-^)d^  =|1f||. 

By  "llm"  we  mean  either  the  formal  limit  If  It  Is  non-zero, 
or  If  not,  then  the  lowest  order  term  In  k. 

The   Uj_  part  of  the  Integral  (Eq.  3)  becomes  In  this 
notation: 

(10)  u^    =    I    d^^   ^(e-B1i^)(^  XB)  II     -    ^lr?i/{T\-^   +   2ii)\\ 

+    (fi^-^J4j|cos2-^(Ti-V   2^)  II    +    [{i'Bx-h^)^^ 
-    (fT^-^)IXB]  II  sin  ^   cos   if{T]^'^   +   2lf)  \\     . 

In  general.  If  V   Is  a  vector,  then  by  Vj_  we  mean 

(11)  V^  =  BX(VXB)  -  V-B(V-B)  . 

For  the  constant  2i/      contribution  we  note  that  If  V   is 
independent  of   I,   then  (cf.  Appendix  D) 

(12)  <^-V  ^>  =  V. 


The  non-constant  t]~        contributions  are  all  of  the  form, 
(V  =  V,) 

(13)       L  =  <(^.V)l(k^.^j2^-2> 


-  55  - 


This  Integral  Is  calculated  In  Appendix  D,  and  is  there  shown 
to  be 

(14)  L(V)  =  2^Ni  k^^"^  L(V)  . 

The  linear  operator  L  rotates  and  stretches  V 
according  to 

(15)  L(V)  =  b2^k(k-V)  +  bgj^  kX  B(k><  B-V)  . 

The  unit  vector  k  Is  in  the  direction  of  k.^,      and  the 
bp^  coefficients  are  given  by: 

(16)  b,^  =  l->5--'{2B-l)       .    ,'   .  b2^(2N-l)-^  . 

2N      2-4-.-(2N)  2N     2N 

There  are  four  distinct  contributions  to  Uj_   of  this 
form.   They  arise  from  the  terms,   (^ -B  x  Ti-^  )^  X  B,  ^ -fi-^^  , 
{i'BX^^)i,      and  (i'-fi^)^  X  B,       (noting  that  ^  =  ej_=  B  x  (|xB) 
only,  and  similarly  for  1i-^ ) .   These  terms,  respectively, 
integrate  to,   L(BX^^)xB,  L('fi^),  L(Bx-fi^),  and  L(1i^)XB. 

Let  us  rewrite  equation  (lO)  as 

(17)  u,  =  u(°)  +  .p)  , 

Where  u[°^  is  the  contribution  from  the  2f  integrals,  and 
u[  ''  the  contribution  from  the  r\~  integrals.  Combining 
the  above  yields 

(18)  u[°)  =  pi [ -IN -f^-,  +  B  X  ^-1  ]  , 

2(N  -  1) 
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^N-2  ,  2N-2  ^,, 

(19)       u[^^    = !  -B  v^  L(B  K  li,  )  +  L(ft,  ) 

nj^(2N-2)l     ;  ^       ^ 


+  1[L(B  Kfi^)    +   B  X  L(-1i^)] 


Prom  the  latter  equation  we  see  that  In  order  to  maintain 
a  finite  phenomenon,  t]        must  be  of  the  form 

(20)      n^   =  k2N-2  [a^+  ka;+  ...]  . 


The  term  o        Is  Independent  of  k  and  Is  to  be  determined. 
To  obtain  a  closed  system  of  equations,  we  must  Include 
the  lon-veloclty  equation.   Since  the  value  of   la   Is  not 
Integral  about   a  =   IN,   a  more  direct  means  of  Integration 
may  be  employed,  which  shall  be  Introduced  in  the  succeeding 
section,  the  results  of  which  we  pre-quote  now  (cf.  Eq. 
C-III:  3^) 


(21)       0^  =  ^^2^— '1^1^  B  +  a^^V  . 
a  +  1  (- 


We  recall  that  cu  is  such  that  a  =  IN  +  t]-^,      so  that  about 


this  frequency 


(22) 


n&  =  (v-v)  +  im  +  o(k2^"^) 


(v-v)/n  -  IN  Vm  +  0(k2N-2) 


An  immediate  eigenfrequency  is  obtained  in  the  said 
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limit  k  — >  0,   If  In  addition  we  let  ^   — >  0,   and  consider 
the  case  of  small  self-collision  frequency  (v),  but  large 
Vj_  (cross-collision  frequency).   In  this  limit,  (cf.  Eq.  C-I; 
32),  we  obtain 


(23)      "^1  =  ^1  =  -  ^3^" 


where   v^  ==  v^  (-^  =  O),   and  again  cd  is  given  by 
a  =  IN  +  ri^. 

Equation  (21)  for  the  ion  velocity  is  in  the  form 
Uj_  =  F  Au,   and  equation  (lO)  for  the  electron  velocity  is 
in  the  form  Uj_  =  PAu,   The  F's  are  the  implied  linear 
operators.    Subtraction  of  these  two  equations  yields  the 
closed  set   Au  =  (F-P)Au.   The  condition  for  a  non-trivial 
solution  yields  the  desired  roots,  viz.,  det|F-F-l|  =  0. 
This  condition  appears  as 


(24)       {z-lf  +   w^ 


where 


(25)      z  =  V   [|-  -  INy  -  Sa] 


(26)  w  =  v^  [1  g-  +  y  +  a]  . 

The  coefficients,  p,  y,   and  a,   are 

(27)  S  =  (a2  +  i)-\   p   .  flhlLJi   ,      y  =  1  (n2-  i)-1 

^     (2N-1)J   2N         2 
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The  solutions  to  the  above  secular  equation  appear  as 
(28)      o^-^    =   0 


and 


(+) 


(29)      ^r'  -  2p^T 

where 

It  follows  that  there  are  two  roots  about  each   a  =  IN, 
and  these  are  given  by 

(31)      Ico^^^  -  V  -  ±m   -   2pNnT  k^^"2 


and 

(32)      loD^")  =  V  -  INfi  -  a^k^^'    , 

where   N'  >  2N-2 . * 

These  relations  Indicate  that  the  first  order  decay 
time  Is  as  the  Inverse  self -collision  frequency,  while  to 
terms  of  higher  order  (i.e.,  terms  of   0(k    )),  the  decay 
time  (i.e.,  [v  -  2^12  Re  t]~  ),  is  a  function  of   a  and 
V-,   through  the  parameter  x.   We  immediately  note  that  if 


That  this  is  so,  becomes  clear  if  we  write  ok 

k2N-2(~  j^N"^^   N"  >  0,   in  which  case,  to  order  k^^"^, 

„   ~   N" 
the  "coefficient"   a  k    would  be  zero. 
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V  =  V-,  =  0   (no  dissipation),  then  Re  x  =  0,   and  there  Is 

no  decay  of  these  magnetic  oscillations.   The  most  significant 

part  of  the  dispersion  Is  In  terms  of  the  frequency   v^, 

which  Includes  the  effect  of  the  number  density  {^     =   en  /e  B  ) 

and  the  cross  collision  frequency  v  .   We  note  that  the 

explicit  orientation  of  c|d   In  the  formulation  Is  such  that 

the  Influence  of  the  density  of  particles  on  the  "pure" 

resonance   Nfi  can  be  eliminated  by  a  large  enough  magnetic 

field,  and  Inversely  the  gross  ordering  effect  of  the  steady 

magnetic  field  can  be  weakened  by  a  large  enough  number 

density.   Since  the  ratio  of  plasma  frequency  to  Larmor 

o 
frequency  appears  in  4)  ,   and  is  usually  quite  large,  the 

latter  of  the  two  situations  would  appear  to  prevail.   It 

should  be  remembered,  however,  that  these  effects  are  higher 

order  in  the  expansion  parameter  k,   than  are  the  "pure" 

resonances.   The  presence  of  the  other  gas  (i.e.,  the  ion 

gas  in  this  case)  is  most  directly  exhibited  by  the  presence 

of  the  coefficient   a,   which  is  evaluated  about  the  electron 

resonance,  so  that 

(35)     S  =  S  (|^  -  i^)  +  o(k2N-a)  . 

An  interesting  result  is  obtained  if  we  note  the  form  that   a 
assumes  in  the  inverse  calculation  about  the  ion  resonances. 
In  this  case,   a   is  calculated  with  co  such  that   a  -   IN  +  t)  , 
and  appears  as 

m  a  =  m  (5zl  _  ill)  +  Oik^N-S)  . 


By  comparison  of  these  two  latter  forms  It  is  easily  con- 
cluded that  the  effect  of  the  ions  on  the  electron  resonances 
is  greater  than  the  effect  of  the  electrons  on  the  ion 
resonances  by  a  factor  which  is  approximately  the  mass  ratio. 

Finally,  we  note  that  these  resonances  are  seen  to  be 
doublets,  with  a  splitting  which  is  of  order   k     .   This 
splitting  occurs  about  each  resonant  frequency,  two  about 
the  ion  frequencies  and  two  about  the  electron  frequencies. 


C-III.   Formulation  of  Generalized  Macroscopic  Equations. 

The  standard  technique  of  forming  the  macroscopic  con- 
servation equations   by  way  of  moment  operation  converts  a 
closed  system  of  non-linear  integral -differential  (Boltzmann) 
equations  into  a  non-closed  system  of  linear  differential 
(moment)  equations.   Introduction  of  the  transport  coef- 
ficients serves  to  close  the  latter  system  and  in  a  phenomen- 
ological  manner  recaptures  the  effect  of  the  collisions. 

The  technique  of  introducing  a  simplified  collision  form 
keeps  the  effect  of  the  collisions  very  explicit  through  the 
presence  of  the  collision  frequencies.   Although  postulating 
these  collision  frequencies  is,  of  course,  ad  hoc,  it  is  a 
statement  in  the  microscopic  realm,  so  that  the  subsequent 


C.  S.  Gardner  has  considered  this  problem  for  the  case  of 
a  collisionless  plasma  and  demonstrates  more  formally  the 
existence  of  such  double  roots,  (report  to  appear). 


61 


analysis  starts  from  microscopic  origins  and  Is  self  contained. 

As  has  already  been  noted,  the  system  of  equations 
(C-II:  l)  Is  a  closed  set  of  ten  homogeneous  algebraic  equa- 
tions In  (k^cb)  space.   We  will  now  embark  on  the  formulation 
of  a  more  explicit  form  of  these  equations.   Since  these 
operations  are  on  the  transformed  variables,  the  related 
coefficients  will  be  functions  of  k  and  co.   Were  we  to 
transform  back  to  (x,t)  space,  a  system  of  Integral  equations 
would  be  obtained  which  would  resemble  the  Integral  Inverse 
of  the  differential  macroscopic  equations  written  In  Section 
C-I,   the  distinction  being  that  the  former  set  would  be  a 
closed  system  with  ten  dependent  variables.   However,  as  we 
shall  see,  a  differential  system  may  also  be  obtained  through 
forming  proper  linear  combinations  of  the  set  (C-II:  l). 

Still  another  choice  Is  available  If  we  view  the  forma- 
tion of  the  transform  equations  In  the  equally  legitimate 
sense  of  seeking  wave-like  solutions  of  the  form 
exp^lk-x  -  lojtl  (i.e.,  as  opposed  to  considering  the  transform 
operations  as  a  mathematical  artifice  which  serves  only  to 
determine  the  solution  In  (x,t)  space).   This  Interpretation 
of  seeking  plane  wave  solutions  permits  the  formulation  of 
frequency  dependent  transport  coefficients. 

Once  again  we  consider  the  equations  (C-II:  l)  and,  con- 
sistent with  the  Investigation  of  long  wavelength  disturbances ^ 
set 


The  pseudo  Inhomogeneous  terms  are  taken  to  be  those  arising 
from  the  collision  forms. 
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(1)   e^ta^  +iy  Bin  ^]  ^  ^25^|^^  ^^^  ^^^  ^  _  2y2^.^2,^  j,  ^  ^(^2^ 


where, as  in  Section  C-II,  we  recall  that  y  =  kj_'^j_'  ^nd  that 
a  =  a  +  Ix,   a  =  V  -  Ico,   x  =  k  -^^^      (note^  x  +  y  =  k-|). 

This  present  analysis  will  include  terms  to  0(k  ). 

Substitution  of  the  above  expansion  into  Eq.  (C-II:  l) 
yields  for  the  first  three  moments 

(2)   I  =  f  dh      i^J   h   d-^  e2^'^[l+  21y  sin  i>  -2y^s±n^1f]^{R^^) 


O)      I  T  =  j  dh      ^J   ^   d^  €^^"^[1  +  2iy  sin  ^ 

2 
-  2y2sin2^]^(R^U(^^) 


(4)   I  =  j  dh      ^g-   A)  d^  e2^'^[l  +  2iy  sin  "^ 
-  2y^sin2^]'h(R^4)R_^e  , 


(with  an  image  set  for  the  ion  moments),  where  we  recall  that 

(5)  'fi(R^I)  =  fiQ  +  fii{B(B-^)  +  (sin^)^xB+  (cos  t)Bx(exB)|, 

= -fi  +^l^„    -    (sin  T/')^-fi-|_  X  B  +  (cos  V')^-fiJ 

(6)  R_T/,^  =  B(B-I)  -  (sin  ^)e  X  B  +  (cos  T/')Bx(^  XB)  . 
Recognizing  that  all  of  the  dif     integrations  are  of  a 
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similar  form,  we  define  the  following  functions  of  a  =  v  -  1cd+  Ix 

(7)       s^C^  =   -^1  6   e^^^  sm^^  cos  N^  d^. 
e   -1  -' 


To  perform  these  Integrations,  we  note  the  Identity 

N+M 

(8)       sin  1/   cos  i/  =  ^       d  cos  ni/  +   b  sin  ni/    , 
n=0   ^ 

and  then  write   sin  nil     and   cos  nif      In  their  equivalent 
exponential  form  to  obtain 

hira.  -,x-l  X     2B.il   „„       ,    ^,,,       r  2a 


(9)  (e^^^-1)"^  4)  e'^^'^cos   nii   6.ii  =  C 


^        ka^  4-   n^ 


(10)  (e^^^-1)-^  j)  e^^^sm  n^  d^  =    ^      " 


^        ^a^  +   n^ 


so  that,  for  Instance, 

A^  =  0  ,      XC^  =  1/25 


(11) 


nC^+   2a^^  =  0. 


A  few  of  the  pertinent  Integrals  are  listed  below: 
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(12) 


If  these  results  are  substituted  Into  equations  (2-4), 
there  remains  only  the  d^i   Integration.   Were  we  to  perform 
these  Integrations  afortlorl  with  the  S  C   function  of  x 
left  Intact,  the  resulting  equations  would  be  valid  to  terms 
of   0(kj_)   for  the  perpendicular  part  of  the  wave  number; 
however,  one  would  obtain  a  closed  formulation  for  the 
parallel  component   k  .   This  asymmetry,  which  Is  Inherent  In 
the  problem,  has  already  been  exhibited  In  the  previous  section 
on  magnetic  singularities  wherein  the  formulation  included 
terms  to  0{k^J^~^),      but  only  to  0(l)  in  k„ . 

An  example  of  the  opposite  extreme  is  obtained  by 
setting  y  =:  0,  and  examining  the  resulting  formulation 
with  all  functions  of  x  left  In  closed  form.  The  classical 
Landau  damping  is  Included  as  a  special  case  of  this  class  if 
we  further  let  all  of  the  disslpatlve  terms  vanish  (i.e.,  all 
collision  frequencies  equal  zero).  In  this  extreme  limit  the 
master  matrix  becomes  diagonal  so  that  the  number  density  and 
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parallel  component  of  the  electric  field  are  Independent  of 
the  other  variables  and  are  governed  by: 


-eE 


(13) 

+eE 

mC 


dhA 


-Ico+TkTC 


dh  0 


■lud+  ik  4    C 


(note:      k     Is   dimensional) 
Polsson's   equation   is 


(14)  ik  E     =  — ^   (n-n)    , 


and  serves  to  close  the  system.   We  recognize  these  to  be 
the  generating  equations  for  the  Landau  phenomena.  '    '     ' 

The  dissipation  associated  with  this  collisionless 
damping  is  understood  to  be  a  diffusion  of  the  charge  accumulc 
tlon  regions  of  the  wave.   We  will  be  concerned  with  the  more 
direct  dissipation  due  to  collisions. 

Returning  to  the  more  general  formulation  we  see  that 

the  S  C    functions  are  linear  combinations  of  ^   and  /L    . 

These  are  all  of  a  similar  form  and  include  the  factor 

[(2a  +  in) (2a  -  IN)]"  .   Dlmensionally  these  factors  appear 

as   [2v-i(2cD  t  m)   +   2ik__|J"^.   It  follows  that  the   d|__ 

part  of  the  remaining  Integration  is  well  defined  if 
*  —  ■ 

Cf.  footnote  p.  (42);  also  to  obtain  the  more  standard 

Laplace  form,  replace  -ico  by   s. 


*-)f 


In  this  manner  the  magnetic  res 'nances  are  kept  in  the  for- 
mulation.  They  become  explicit  after  the  d^,,  integrations 

are  performed,  and  appear  in  the  [S  C  ]'^  coefficients. 
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Re<|2v  -  21cu  t  INfij  f   0.   This  condition  Is  assured,  for 
Instance,  If  we  restrict  the  Investigation  to  roots  in  the 
domain  Re(-lQ)  >  -v.   Such  roots  describe  phenomena  which 
decay  no  faster  than   e""^  .   The  expansions  now  presented 
are  asymptotic  to  these   d^    integrations. 

If  the  d"^  integrals  in  equations  (2,3,^)  are  replaced 
by  the   S^c'^  functions  (12),  there  results 

(15)    g  =^c°|fi^+ihje„l  -sji-fi^KB]  +  c[^hJj+  s[2iy(fiQ+fi;;^jj 

-  s2/21y(e-fi^xB)|    +   Sc[2iy   ^fij|    -   ^^[2/ {^ ^+ "^[^ ^;)\ 
+   S^{2y2^fi3_KB]   -   S2c[2y2|.ft^|      =    (V^ 

(16)  Jt=<V(^)) 

(17)  ^  =<;^c°|(fi^+fi"^J4„J     +    s|^j2iy(fi^+'fi;;iJ-|-h^XB] 

+   Sck_(21yMT^)    -    ^XB(^ti^)    +    ^J2ly(li^  +  1i^lj 

-  e-ti^XB]]   +    s2|-^j2y2(TT^  +  ti;;^J  +  2iy    l-li^XB] 


More  generally,  the  integrations  are  well  defined  if  either 
V  +  IM(co)  ^   0,   or   2Re(Go)  +  Nfi  +  2x  =?s  0. 
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+   C^||^4-1T^/+   S^^„(2y^^Ti-^><B)    +   ^X  B[2y2  (TT^ +1ii;e  J 
+   21y   l-Ti^XB]]  +   SC^|e^21y   e-?i^]  +   CS^U  J -2y^  ^ -Ti^) 


-   |^[2y2(Ti^+1i^^J   +   21y   ^Ti^KB]    -    ^  ><B(e  •1!^21y ) 
+   S^f-^XB   2y2Mi^XB]+   cV|4^^!T^(-2y2)| 
+   S^c{l^(2y2Mi3_XB)    +   ^XB(2y2^1^J)|^         . 


Many  of  these  Integrals  are  annihilated  immediately  by 
odd -even  arguments.   The  remaining  terms  appear  as: 

(18)  I  =(c°[[^o  +  -Rie„]^„|  +  Sc[^_^2iy  M!^-^XB(M=i^) 

+  ^J2iy(li^+ti^^j  -  ^tl^XB]]  +  S^l  -2B[^^y^^^ 

+  ^_^iy  ^tr^XB]  +  ^XB[4-h^XB  -  2iy (TT^  +  1^^^  J  ]] 

+  c2||^e.-fi'|  -  S^^I^XB  2y2M^^XB]-  c2s2[2y2|_,  ^TT^  ] 

+  S^C  •^^2y^^Ti^XB  +  4XB  2y^^1^^  |  \ 

(19)  I  =^C°|tT^  +  -R![^„]-  s2|2iy  Mi^Xb]  +  Sc[21y  ^tl^^  | 

-l^^^-ol) 

(20)  ^  =(^V^  tC°K+'^l^„]  -  S^l^^y  ^-H^XB  ] 


SC  [2iy  4-11^1  -  S2^2y-  ti^]]^ 
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The  S  C   functions  of  x  are  expanded  as 


(21)    sV=  [cV](°)-l4sV](l'-x2[sV](2)  H.  ...  , 


Which  serves  to  define  the  coefficients   [S^c'^]^'^^. 

Recalling  the  definition  of  the  S^c'^  functions  (cf. 
Eq.  7),  we  see  that  the  relevant  coefficients  involve  the 
following  four  forms: 


(i^  i)-l=(a2+i)-l-ix^^-x2X5a!zll^  . 
(a2+l)2      (a2+l)^ 

_  ix  (^!-  n    -   x2  a(a^-3)  ^  . . . 


a^+   1  a^+  1      {a.^+lf  (l+a^)^ 

(22)  ,   2    ^       -  ^ 

1  ^    1     _  ix   (3a  +  1)  _  ^2  ba  +  3a  +  1   ^ 

a(i2+l)   a(a2+l)      a^Ca^+l)^       a^(a2+  l)^ 


(2a)~^  =  (2a)"^  -  ix(2a2)"^  -  x2(2a^) 


The  related  coefficients  appear  as 


[C°](0)M2a)-\   [s2](°).[4a(a2.l)]-\  [c^  ]  (-)  =  ^^2^) 

[SC]^°^  =  -[4(a2+l)]-\   [C^-  s2](°)  =  a[2(a2+l)]-\ 

(23) 

[s5c](°)  =  -g  [(a2+4)(a2+l)]-\   [s2c^']^''''=[l6a(a^  1 )  ]  "^ 


[s'^]^^)  -   3[^a(a2+l)(a2+i+)]-\   [SCj^^^  =    '%  o 

2(l+a  ) 
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(23) 

ts2](l)._^^,   [C°](2).(2a^)-\   [c2-s2](2).a(:a!^ 

2.(2)  ...  6aV  3a^+  1 


:s^] 


^a^(a^+  1)^ 


qpl(2)  _  _  1  I    3a^+  1 


Ill-a.   The  energy  equations. 


Substituting  these  forms  into  the  above  equations,  the 
<  >  integrals  (i.e.,  d^^ )  are  readily  performed  (cf.  Appendix 
D).   Let  us  first  consider  the   n  and  T  equations: 

(24)  I  =.  vn[C°](°)  -  k^[C°](2)(^^  ^^)  ^  2l[s2](°^  kxB 

+  2i[SC](°)k^-fi-^  -  2k2[s2](°)(vn  +7)  -i[C°](^)k_^'h^ 

(25)  f  =  |7[C°](°)  -k2[C°](2)(vn  +  Ir) 

+  2l[s2](°)fi^-kXB  +  2l[SC]^°^k^-1i^ 
-  2k2[s2](°)(vn  +  i7)    -   2l[C°](l)k  hi' 
where  we  have  set 


y   -  vT  +  VgAT  =  (v-V2)T  +  v^T 
(26)    ^    _ 

y  =  vT  -  VpAT  =  (v-v„)T  +  v^T 


Cf.  footnote,  p.  (HO). 
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Substraction  of  the  temperature  equation  from  the  number 
density  equation  yields,  (returning  to  dimensional  form) 

^/^_5n)_„   AT  ,  5  .2/  ^^         ,  ^n    "^ 


o 


(27)    ^  ^  -4^  =  V,  ^+2c-  ^^^+_1^^^. 


^t   P^    3  n^     -'2  T^  ^  3  "     2^„2  '   2   1   T 


o   ^    I  a^+$7^   a^   j   ^o 


rP7  n   ^  1  P     5  n  1    ^^   AT  _^  3  ^2  I   "^^    ,  ^"    5! 


Again  we  have  partially  Inverted  back  to  (x,t)  space,  which 
enables  us  to  obtain  In  this  Instance  an  00-dependent  heat 
conductivity. 

If  this  equation  is  compared  with  the  ellmlnant  of  the 
Individual  energy  equation  (B-II:  16)  and  the  continuity 
equation  (B-II:  j) ,    viz., 


(28)      §t 


^,  P_  _  5  n_\  ^  V-C    ^    ^      AT 


\Po    5% 


2  T 


o 


one  obtains  the  following  expressions  for  the  components  of 
the  Individual  electron  heat  conductivities: 

r-      w  C^   _,   I  VT+  v^At)    ^ 
■^    -^  a  +  n  I     o    J   m+m 

(29)  2   r 

^  w^C       VT+  V  AT  /    ^ 
"     -^  a    "      o       i^+ni 


Similarly  for  the  Individual  Ion  heat  conductivities  we 
find 
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(30) 


^  a  +  n  T„        m+m 


^  w^C^      vT  -  v^AT  } 


T     1    m+m 


The  more  relevant  total  heat  flow  is  merely  the  addition  of 
the  above  forms  and  is  given  by: 


-^   o      a  a  J      \m+m/ 


m+m. 


(31) 


f  ^  \7x^  ^|-^(^T  -  v^AT)  +  -2^(vT  +  v^AT) 


o 


a  +f^  a  +f2 


"^-"^  >  w  Au, 
o   -» 

m.+m 


In  the    limit    of     m  =  m,      T  =  T     and      v  »  03,    the   above   yields 

f  1  8 ") 
the   Chapman-Cowling   result, 

(52)       ^  =  7—72       • 

^M       1+   (f^/v)"^ 

More  generally,  if  AT  ^   0,   we  see  that  the  heat  flow 
is  proportional  to  the  gradient  of  a  linear  combination  of 
both  temperatures.   In  addition  the  contributions  from  the 
difference  in  flow  terms  indicates  that,  even  to  the  lowest 
order,  the  gas  mixture  appears  non-adiabatic .    ' 

-  72  - 


Ill-b.   The  continuity  equations. 

To  obtain  the  continuity  equations  we  form  a  linear 
combination  of  the  number  density  equation  (24)  and  the 
velocity  equation.   They  appear  as: 

(33)  I  -  vn[C°](°)=:-k2[C°](2)(^^^^^^2i[s2](°\.k><B 

+  2i[SC](°^k.fi^  -2k2[s2](°)(vn+;r)  -l[C°](l)k_;ft^  . 

(34)  ^=  -  ikJC°](l)(vn+r)  +  fti[C°](°^  -3A;[C°](2) 

+  2[SC]^^^k_^k^-h^  -  2h^xrB[SC](°)  +  2h-^><;  B  k^[sc]^^) 
+  2ik[SC]^°)(vn  +  ^)  +  2h^[SC]^^)k_^k,-  2k2f,^[  S^  ]  ^°) 
+  2[s2](l)k^^fi^-kXB  +  fi^[c2-s2](°)-  k2fi^[c2-s2](2) 

-  2[S^]^°^fk5ft^  +  2(k  XB>r-i^-k  X  B? 

-  2[C^S^]^°)  |k^h^  +  2k^k^-h^j+  2[S^C]^°^|-2k^BXh-^ 

-  2kj_h^-k  X  B  +  2k  X  B  h^-kj_ 


The  linear  combination  is   l/2(vn  +  ik-u),   and  to  0(k  )  Is 
(35)    I  (an+  ik-u)  -  vna[C°] ^^^ =  k2(vn+ ^  )/[C°] ^^ ^ -a[C°] ^^H 
-  lk,,ti;;|  -a[C°]^^^  +  IC^](°)]  -  2[ik.-ri^X  B 
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+  k2(vn  +  r)l|a[s2](°)  +  sc(°'|  +  Ik-h^  (  2a[SC]  (°) 


c2-s2](°)   +  0(k5) 


Reference  to  the  table  of  functions  on  p.  69  ,  Indicates  that 
all  of  the  terms  on  the  right  hand  side  vanish  identically,  so 

o 

that  to  0(k  )  we  obtain  the  standard  continuity  equations: 


con  -  n  k-u  =  0 


(36) 


It  should  be  noted  in  passing  that  although  these  equations 
are  here  obtained  to   0(k  ),  they  were  also  derived  in  Section 
B-II  through  the  more  standard  procedure  of  moment  construc- 
tions, from  which  it  follows  that  they  are  valid  to  all  orders 

in  k. 


III-c.   The  momentum  equations. 

To  Invert  the  velocity  equations  we  change  from  a  vector 

representation  to  a  cartesian  representation  which  is  defined 

with  respect  to  the  left-hand  triplet  of  unit  vectors 

[kj_,  k_LXB,  B].   The  system  is  such  that  if  V  is  a  vector 

then  V^  =  B-V,   V^  =  k-V,   V  =  kXB-V.   In  forming  the 
^  X  y 

three  components  of  the  velocity  equation  we  note  the  following 
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identities,    (calling     kj^,    k) 


k-h-j_x  B  =    -h^-kx  B  =   -hj 


(37)  kXB-h-^XB  =  h^-B  >^(kXB)    =  h^-k  =  h^ 

k  KB-h-^   =  h^ 


k-k  =  k^,        kxB-kxB  =  k      ,        k     =0 


The  three  components  of  equation  (C-III:  3^)  appear  as 


(38) 


^=  -ikJC°](l)(vn+7)  +  n,^|^[C°](°)-  3k2[C°](2) 

-  2k2[s2](-)]  +-h-[2k^k^[SC](^)j  +1^j(2k^k^[s2](l)j 

(39)  ^-  2ik^[SC]^°)(vn+^)   fi^[2k^k^[SC]^^)] 

+  fi^[[c2-s2](°)-k2[c2-s2](2).2k2([S^](°) 

+  S[cV](°))j  +-fiy[2[SC](°)-2k2[SC](2)_8k2[S^C](°)j 

(40)  ^-  -21k^[s2](°^vn  +  7  )^--n^[-  2k^kjs2](l) 

+  ^j[-  [SC](°)  +  2k2[SC](2)  ^8k2[s5c](°)| 

^-r.j[[c2-s2f-\2[c2-S-](2)-2k^(3[S^]^°^[cV](°))], 
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If  the  u       equation  Is  multiplied  by  a  and  then  added 
to  the  u   equation,  there  results: 


(^1)      -^u  +  au^  +  :^  (vn  +  ;^  )  -  ^J  +,^^^1  =  0 


The/T^^.  matrix  Is  a  function  of  the  [C^S^]^^^ 
coefficients. 

In  similar  manner  If  the  u   equation  Is  multiplied 
by  -a,  and  then  added  to  the  u  equation,  we  obtain 


(42)      -u^+  auy  +  h^  +^5^y^h^  =  0 


Finally,  the  z  equation  Is  rewritten  In  the  similar 
manner: 

Ik 


(43)      au^  +  _£  (n  +7^)  -  h^  +,/4^/i  =  0 


These  three  equations  are  the  components  of  the  more 
familiar  vector  equation  (return  to  dimensional  notation  by 
multiplying  by  fiC;  recall  that  Q.  =   eB/m) ,  where  we  have 
recalled  that  (x,  y,  z)  is  a  left-hand  triplet  and  that 
k  =  (k^,  0,  k^): 

^^^)        -o  It  +  ^  (E  -  -xs)  -  i^ ^p  -Ay\ 


-  76  - 


similarly,  for  the  ion  equation  we  have 

ri'^    en  ^         a        ^ 

m  am  '^ 

-  -V,  ^  n  Au  +  ^  ^  V(AT)  . 
a  m 


We  have  written  p   for  K(n  T  +  T^n),  cf.  Eq,  (B-I:  5). 

The   a~   operator  can  be  interpreted  in  either  of  two 
ways.   First,  let  us  call   va~  p   the  augmented  pressure   p, 
so  that  (recall  that  these  transformed  variables  have  carried 
an  unwritten  star  --  cf.  p.l42), 


(46)      va-^p^ 


The  inverse  of  this  equation  appears  as 

X 
(47)      p  =  p(0)e-'^^  +  V    e"''(^"^')p(t')dt 

o 


(48)      p  =  P  +  e-^^[p(0)  -  p(0)]  +  J  e-^(^-^')p^,(t-)dt'  . 

o 

The  first  term  is  the  standard  pressure,  while  the  second 

indicates  the  influence   of  the  initial  pressure,  which  is 

seen  to  be  very  small  after  only  a  few  collisions.   The  last 

term  is  a  non-local  effect  over  past  pressure  which,  for 

large  collision  frequencies,  becomes  almost  local. 


It  can  be  seen  immediately  from  Eq.  (46)  that  as  v  — >  co  , 
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The  second  interpretation  comes  about  when  we  note  the 
identity 

(49)    v(v-ia3)-^  ^ll  +  (ia)/v)  r    ...    +(icD/v)"-^|  +  "^^^{^^"   • 

This  expression  serves  to  change  the  integral  operator  to 
a  series  of  differential  operations  with  a  high  order  integro- 
dlfferentlal  remainder. 

In  the  limit  of  finite  wave  speeds  (MHD),  to  be  discussed 
momentarily,  (ci^/k)  — >  const.   Since  the  gradient  operator 
already  multiplies  the  augmented  pressure,  we  see  that  in  this 
limit 


(50)      Vp  =  Vp  +  O(k^)  . 


It  follows  that  in  first  order  theory  the  pressure  term 
becomes  the  standard  kinetic  pressure  term  with  p=K(nT  +  Tn  ) 

Referring  to  equation  (27),  we  see  that   AT  =  C(k). 
Since  there  is  a  gradient  operating  on  this  inhoraogeneous 
term  in  our  momentum  equations,  it  follows  that  the  related 
term  is  of  O(k^) . 

We  note  that  to  0(k),   h-,   has  been  eliminated  in  the 
above  equations  (44,45)-   This  was  done  in  order  to  recapture 
the  more  familiar  form  of  the  momentum  equations,  although  by 
the  analysis  of  Section  B-II,  we  already  knew  that  these 
equations  were  contained  in  the  model.   Since  the  forms  with 
h-^   kept  explicit  and  those  with  h-,   eliminated  (in  favor  of 
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E,  u,   u)  are  related  by  linear  combinations,  the  resultant 
dispersion  relations  of  the  two  schemes  are  identical.   With 
this  point  in  mind  the  above  set  (27;  36;  4^,45)  is  taken  as 
our  final  form  of  the  velocity  moment  equations  to  0(k  ). 

We  now  list  the  elements  of  /X)  .  .  to  note  the  explicit 
dependence  on  cos ^  =  B-k/Bk  and  for  reference  in  higher 
order  calculations.   In  addition,  if  we  compare  the  momentum 
equations  of  this  section  with  equations  (B-II:  12,13) j  we 

see  that  the  elements  of  yy  ^ .  enter  inthe  identification 

■x- 
of  the  stress  and  electrical  conductivity. 


(51) 


V.  L 


l{<i  -  <X.} 


1  r-^"-^ 


2k  k  ^ 
X  z 


-'J^< 


The  orientation  is  in  the  order  x,y,z,  with  B   in  the 
z  direction  and  k  in  the  x  direction.   The  submatrices  are 
given  by 


(52)    ^ 


^{2[SC]^2)-a[c2-s2](2)j     .{2a[SC](^).[c2-S^](2)) 
i-{2a[SC](2)^  [C2-S2](2)}    .(2[SC](2)-a[ 


(2)  .rc2_s2l(2) 


^  Recall  that   hi   is  composed  of  three  parts.   The  E  and  Au 
terms  contribute  to  the  conductivity,  while  the  remaining  u 
(and  Au)  terms  contribute  to  the  viscosity. 
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2 


2a[S 


(o)} 


\ 


[8a[s'o](°'  +  2[S^](°)      -/8[s'o](°'-  eats"]'"' 

+  6[sV](°)|  -  2a[cV](°)]  / 


2 


Aa[scf+   [S2](l)) 


SC](1'+  a[s2](l) 


P  =  ■••  {3t0°] 


0,(2) 


^=({[so](i)|   ([S^] 


,    g'=..{2[S-](°)i  . 


The  explicit  dependence  on  a   of  all  of  these  coefficient! 
appears  on  p.  69  . 

The  aggregate  of  equations  which  we  have  called  the 
energy,  continuity,  and  momentum  equations  are  seen  to  form 
a  closed  set  of  linear  equations  in  the  ten  unknowns  (n,  u,  T; 

/N     /N     /S  .  /   2  \ 

n,  u,  T).   These  equations  are  formulated  to  0(k  );  however, 
their  co-dependence  is  in  closed  form. 

In  the  succeeding  section  we  will  consider  the  various 
roots  of  the  related  Master  determinant.   We  shall  see  how  the 
spectrum  of  plasma  phenomena  is  generated  by  examining  the 
different  forms  of  solutions,  cu  =  ^^(k). 
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C-IV.   Dispersion  Relations. 

It  Is  convenient  at  this  point  to  list  the  system  of 
equations  derived  In  the  previous  section  which,  together 
with  Maxwell's  equations,  form  a  closed  aggregate: 

(2)  ||+n^VG   =   0 

*"  ^\Po       3"o)^     2T„+  3°    ja^+Q^       a^     I    ^o 

,4,  S     J'P  SS    )_      ^      AT  ^   5  g2i_Zi_  +  Zi]     2^ 


(5)  "ol^^^e^  +  ^x^)  ■^SsVP+Aj-'^l 


(6)  n^  f   -  ^   (E  +   GXBJ    +  ^Vp  +/f^ij-h-L 


=    -V      -     n^Au  +  -^^    \7(AT) 
^  m  am 


(7)  VXB  =   H^J  =  E/c^ 


(8)  \7XE  =    -B 


(9)  V-E  =  Q/^o 
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(10)       V'B  =  0  . 

The  presence  of  the   a  terms  symbolically  represent 
integral  operators  over  past  time.    The  essence  of  these 
operators  Is  such  that  for  large  collision  frequencies  (i.e., 
large  compared  to  the  disturbance  frequency),  they  become 
extremely  local.   This  is  seen  most  directly  from  the 
coefficients  themselves,  since  a  =  v  -  ico.   An  equivalent 
statement  which  follows  from  examining  the  related  integral 
inversion  is  that  the  terms  become  local  after  a  sufficient 
number  of  collision  times. 

Let  us  consider  the  order  in  cd  of  the  transform  of  the 
above  system  in  the  limit  (v,  v  »  co) ,  so  that  all  coefficientE 
are  local  in  time.   This  set  is  seen  to  be  a  system  of  equa- 
tions for  the  variables  (n,  n,  u,  u,  T,  T,  B,  E).   Each  of 
the  first  eight  equations  is  a  first  order  differential 
equation  in  time  for'  each  of  the  above  eight  variables, 
respectively.   The  remaining  two  equations  {9,    10)  serve  to 
eliminate  one  component  of  E  and  B  respectively  in  terms 
of  other  variables  of  the  system,  so  that  only  two  of  the 
three  components  of  the  vector  equation  (7)  and  two  of  the 
three  components  of  the  vector  equation  (8)  are  relevant. 
The  remaining  two  components  are  redundant.   If  we  envision 
the  related  coefficient  determinant  of  the  transformed  system. 


Symbolically,  in  the  sense  that  the  presence  of  the  inver- 
sion integral  is  tacitly  assumed. 
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the  fourteen  diagonal  elements  each  contain  an  to  entry. 

It  follows  that  the  related  dispersion  relation  is  1^^^  order 

25 
in  o).   Oster  -"   considers  a  system  of  equations  similar  to 

the  above  set,  except  that  they  are  temperature  independent, 

and  do  not  include  the  effect  of  collisions.   The  Maxwell's 

equations  part  of  this  system  excludes  the  two  divergence 

equations.   The  related  determinant  appears  l6   order  in  cjo. 

But  we  have  already  seen  that  such  a  system  should  yield 
a  \K        order  relation.   In  addition,  the  related  system  implies 
that  the  time  derivatives  of  the  divergence  equations  are, 
respectively,  zero.   Since  the  divergence  equations  are  not 
included  in  this  system,  two  of  the  fourteen  roots  must  be 
zero,  so  that,  in  all,  four  of  Oster 's  sixteen  roots  should 
degenerate  to  zero. 

Returning  to  the  above  system  (l-lO),  if  terms  of  0(k:) 
are  retained,  but  the  restrictions  on  the  self-collision 
frequencies  (v,v)  are  removed,  then  specifically,  the  trans- 
formed momentum  equations  now  include  an  cd  factor  which 
multiplies   u,   and  u,   respectively.   It  follows  that  in 
this  more  general  case   the  dispersion  relation  is  l6th  order 
in  CD. 

The  higher  order  /f)  .  .   and  AT  terms  of  equations  (3-6) 
can  be  seen  to  contribute  to  the  stress  by  comparison  with 
equations  (B-II:  12,13).   However,  these  terms  are  not  in 
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standard  form  ,  and  are  more  useful  In  the  Investigation  of 
higher  order  dispersion  relations. 

We  will  concern  ourselves  in  this  section  with  the 
lowest  order,  non-trivial  dispersion  relations  which  emerge 
from  these  equations. 

By  transforming  back  to  (kjO))  space  an  algebraic,  homo- 
geneous set  of  ten  equations  in  the  said  ten  unknowns  is 
obtained.   The  condition  that  a  non-trivial  solution  exists 
(viz.,  that  the  master  determinant  vanish)  yields  the  desired 
relations,  which  may  be  written,  03  =  co(k). 

Consistent  with  the  expansions  involved  in  deriving  the 
above  set  of  equations,  we  similarly  expand  the  dispersion 
relation  about  small   k  (cf.  Appendix  B), 


(11)       03  =  00^°)  +  ka3^^^  +  k^J^^ 


(k  is  non-dimensional). 

The  speed  of  light  is  Included  in  the  following  manner. 
It  is  either  taken  to  be  constant  with  respect  to  the  limiting 
process  k  — ^  0,   or  it  is  taken  to  be  0(k)~   in  this  limiting 
process.   In  the  latter  case,  the  parameter  -^  =  ck   is  con- 
stant as   k  — >  0.   Since  this  parameter  ^   enters  in  the 
Maxwell's  equations  part  of  the  system,  which  have  been  kept 
intact.  It  follows  that  the  resulting  formulation  is  closed 
in  the  variable  ^  . 


By  "standard"  we  mean  a  form  from  which  the  viscosity 
coefficient  may  readily  be  identified. 
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The  series  (ll)  yields  two  distinct  classes  of  roots, 

according  to  which  od^°'^  =  0,   or  co^°^  ^   0.   The  distinction 

between  these  classes  is  physically  relevant.   The  first 

categorization  is  indicative  of  finite  wave  speeds 

[lim   (co/k)  =  cu^  '  ^  .      The  latter  is  characteristic  of  phenom- 

k  ->  o 

ena  associated  with  infinite  phase  speeds,  since  in  this 

instance   lim   (oj/k)   is  unbounded.   Let  us  consider  the 

k  -T'  o 
parameter  c   with  respect  to  both  these  limits.   In  the  first 

instance  (finite  wave  speeds)  the  ratio  (cjj/k)/c   is  zero  if 

c  =  0(k~  )   and   finite  if  c   is  constant.   In  the  situation 

wherein   (co/k)   is  unbounded,  the  ratio   (a)/k)/c   is  finite 

in  the  c  =  0(k""^)  {l.e.,yp.   =   const)  limit. 

Returning  to  the  series  (ll),  we  recall  that  in  Section 

C-II,  we  have  already  obtained  one  class  of  roots,  which 

results  from  setting   k-B   =  0.   These  roots  are  of  the  form: 


(12)       .  ^  4°)  .  a„.2N-2 


N       N 


where  uj^°^  is  identified  through  3.{oi^')    =   iN  (cf.  Eq .  C-II:1') 

N 
These  roots  are  such  that  there  is  a  doublet  about  every  o)^, 

so  that  this  class  represents  a  four  fold  Infinity  of  roots. 


Cf.  footnote  p.  50. 
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IV-a.   Magneto -Hydro -Dynamics:   MHD. 

We  will  now  consider  in  more  detail  the  roots  that  are 
Included  In  the  finite  wave  speed  class.  I.e.:   03^°^^  =^  0. 
Referring  to  Appendix  B,  we  recall  that  the  technique  of 
solution  envisions  an  expansion  In  powers  of  k   of  the 
Master  determinant  |A. .|   and  the  system  vector  X^ . 

The  expansion  of  the  vector  X.   explicitly  Involves 


U  =  U^o)  ^  ~y(l)  ^  ~2y(2) 


(13) 


J  =  J 


(o) 


B  = 


:(0) 


E  =  e(°)  + 


f=p 


(o) 


Q  =  Q 


(o) 


(o) 


T  =  t(°^  + 


and  we  recall  that. 


(Note  that  except  for  the   T  terms  we  have  transformed  from 
the  Indlvudual  flow  variables  to  the  "total  fluid" 
variables.)   We  are  now  Interested  In  obtaining  the  lowest 


order  non-trlvlal  form  of  the  system  of  equations  (l-lO). 
To  these  ends  we  perform  the  expansion  indicated  in  Appendix 
B.   Starting  with  Eq.  (7),  we  have 

(7.1)      ik  X[B^°)  +  KB^^K    ...]  =  u  (j(°)  +  kj(l)  +  ...) 


c 


Equating  powers  of   k  yields 

(14) 

IkKB^")  =  H  kj(l)  -  i^E(°)  . 

(We  have  written  co  for  co^  ^  .  ) 

In  similar  manner  we  obtain  the  lowest  order  terms  of 
the  remaining  relevant  Maxwell's  equations.   They  appear  as 


(15)  ik-B^°^  =  0, 

(16)  k  >^E^°^  =  cDkB^°) 


Addition  and  subtraction  of  equations  (l)  and  (2),  respectively; 
yields 

(17)  Q^°^  =  0 

(18)  -kcop^°)  +  Pok-U^°)  =  0. 
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The  addition  of  the  momentum  equations  give 


(19)      -icok   ^U^°)  +  lk-P^°^  =  kJ^^^X  B^ 


and  the  subtraction  gives 

(20)  E^°)  +  U^°^K  B^  =  0. 
Addition  of  equations  (3)  and  (4)  yields 

(21)  pi^,5N(^   , 

o     -^    o 

or  equlvalently , 

(o)  _  10  ^^o  n{o) 


(22)       P^°^  =f  T.P 


To  obtain  the  latter  form  we  have  employed  the  kinetic 
definition  P^  =  N^kT^  (cf.  Eq.  B-I:  ?)  and  Q^°^  =  0  (cf. 
Eq.  17). 

Subtraction  of  the  energy  equations  indicates  that  to 

2  (2) 
terms  of  order   k,   AT  is  zero,  1  .e .  ,  AT  =  k  AT*^  '    (l  +  ...) 

The  related  frequency  is  a  single  mode  and  is  a  function  of 

second  order  variables,  i.e.,  (AT^  ' ,    An^  '). 

The  related  lowest  order  terms  of  the  above  equations 

now  appear  as  (calling  kJ^  ' ,    J;  E^°  ,  E;  etc.). 


Also  note:  2  O  =  m\  +    (m-m) {O/ e) 


(23) 


UXB^ 


B  =  V  X  (U  XB^),     \7*B  =  0  . 


We  have  eliminated  E  through  E  +  U  X  B^  =  0,   and  J 
through  VxB  =   \i^J  +   E/c^,   and   a^  =  —■   kT^M.* 
The  related  dispersion  relation  appears  as. 


2 
(24)       [A'^cos^e  -  v2][v^-(a;^+  A  '  2)v^+a^A '^  (l  -  ^)cos20  ]  =  0, 

c 

We  have  written  v   for  (a/k),  and  A  Is  the  Alfven 
speed,   A   =  ■^o'^^oro'  ^^'^^      ^°^  0  =  B  •  k/B  k.   The  primed 
speeds  are  such  that,  for  Instance, 


c>  A"^ 


The  above  dispersion  relation  Is  referred  to  as  MHD 
(Magneto-hydrodynamics  plus  displacement  current)  In  Fig.  1. 
The  left  factor  Includes  the  generalized  "Alfven"  waves,  while 
the  right  factor  Includes  the  generalized  "fast"  and  "slow" 
waves . 

The  pressureless  MHD*  dispersion  relation  follows  from 


This  Is  the  adlabatlc  sound  speed  for  a  two  component  gas 
with  equal  number  densities  and  temperatures.   With  y 
taken  to  be  5/^,    the  effective  mass  Is  the  average:   M/2 . 
(Cf.  footnote,  p.  38  ) . 

Similar  relations  appear  In  references  (28,  29). 


2 
the  above  In  the  limit,   a  ^>  0.   There  result; 


(25)      v^(A'^cos^e  -  v^)(v^-  a'2)  =  0 


We  have  labeled  the  last  factor,  the  "cold  fast"  wave. 
The  more  familiar  relation  follows  from  Eq.  (24)  In 
the  limit   G  »  [v,a  ,A].   It  appears  as 


(26)       (A^cos^e  -  v^)[v^-  (a^  +  A^)v^  +  a^A^cos^S]  =  0 


1-4 
In  accord  with  the  literature    ,  this  latter  relation 

Is  referred  to  as  MHD  In  Figure  1. 

The  relation  (24)  Is  seen  to  be  6   order  In  (a/k)  . 

Together  with  the  temperature  mode  discussed  above,  these 

relations  Indicate  that  seven  of  the  fourteen  roots  of  the 

system  (l-lO)  behave  as 


(15.1)    CD  -  ko)^^^  +  k^oo^^^ 


The  behavior  (i.e.,  explicit  first  order  terms)  of  the 
remaining  seven  roots  will  be  uncovered  In  the  succeeding 
section. 


IV-b.   Plasma  Electro-Magnetics:   PEM. 

We  proceed  now  to  the  second  class  of  roots  which  is 
characterized  by  the  expansion  o)  =  u:,^^'   +    ...,   cd^°'^  ^   0.   As 
was  already  mentioned,  it  Is  relevant  in  this  limit  to  intro- 
duce the  parameter  -p  =   ck,  together  with  the  constraint 
lim  -f^  =  const. 

We  consider  first  the  ellminant  of  Eqs.  (7)  and  (9)  and 
the   Q,,J   continuity  equation  to  obtain  (cf.  Eq.  C-I:  27), 

(.8)      e„E  =  l^  -"'-/j^^^ 

(calling  CD^  '    =  cx))  . 

The  lowest  order  addition  and  subtraction  of  Eqs.  (5) 
and  (6)  appear  as,  respectively. 


(29)      folt^^JX^c 


and 


(30)       1^  J  -  E^a.2(E  ^  uXB^)  +  n  JXB  =  -v^, 


B   is  a  unit  vector  in  the   B^  direction. 

E  and  U  are  elimated  from  the  (transform  of  the)  last 
three  equations  to  obtain. 


Many  of  the  relations  which  follow,  in  the  absence  of 
collisions  (v,  vi,  V2  ^>-  O),  reduce  to  forms  which  can  be 
found  in  the  literature.   (Cf.  Ref.  10,  20,  25). 
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(51)       bJ  +  1  ^  [BX(JXB)]  +  n   J>CB  -  b'Y^(-^.J)  =  0 


where 


(32) 


2 

l03    03 
b    =    V^     -    Ico   +                °       = 
03     -^^ 

=    03    + 

103    03^ 

b'    = 

-^0              . 

/    2         2x     ' 
03(03   -  ^    ) 

^  -  .^  -  ^'^°" 

n^  = 

=     !'^°2    (m-m)2  ■ 

.   Y?2= 

m  m 

(m  m)"' 

n  =  n  -  fi. 

7  = 

(m-m)2 

(32-) 


which  also  serves  to  define  the  mass  ratio  y. 

The  remaining  continuity  equation  for  _P>U  indicate 
that  in  the  series  f  =  p(°)  +  k  f^^^  +  ...,  0°^  =  0, 
i.e.,  to  lowest  order  the  mass  density  is  not  perturbed. 


Together  with  the  energy  equations,  this  Information  gives 
similarly  that  in  the  related  pressure  series,   P^  '^  =  0. 
The  subtraction  of  the  energy  equations  yields 


(33)       (1^  -  2v^)AT^°) 


This  equation  is  seen  to  be  Independent  of  the  preceeding 
complement.   The  related  mode  is  a  pure  decay  of  the  form 

_ 

Also  called  lower  hybrid  frequency. 
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(34)       At(°)  :=  AT^°)  e"^'2^ 


The  subscript  zero  indicates  an  initial  displacement  from 
the  equilibrium  situation.   In  addition,  reference  to  equa- 
tion (G-III:  25)  indicates  that   T^°^  +  T^^''  =  0. 

Returning  to  the  above  current  equation,  we  now  consider 
the  remaining  3x3  block  of  the  Master  determinant.   To 
evaluate  it  we  choose  a  representation  such  that  ^  =(-^  ,0, -^^) 
B^  =  |B^|(0,  0,  B);  so  that  ^-J  =  fy.J y.   +/'z'^z'   ^^®  secular 
equation  appears  as 


(35)      b(a^  +  Q^)  -  b'[ab^^  ^^f  z^*^  +  "^-  ^^^^  "  °' 


Two   other  convenient   forms   are 


2^     ■    hTf-   b-^2    [i|L+   ^2^    =   0 


(35')  ab(at-b'^    )   +  hn^-  b'^ 

and 


(35")     (a^  +  n2)(b-b'^2)  _  g,bt),(^2_  0=0, 
and  finally,  in  the  form  of  Astrom^^^^  and  Allls^ 


(33,,,)      ,,,2,  ^  (b./^-b)(.^.^^) 
b[ar(a-b'^2)^  ^2^ 


where   0   is  the  angle  between  /O  and  B 
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The  variable  ac     Is  given  by 


(36)      a  =  b  +  1  ^ 

^   ^  CO 


so  that 


1(^2+  03^) 

a-b  '  wD   =  CD  +  


(57) 


b-b 


2    (Ico^  +  oxd) 


Substituting   the   values   of     a:,    h,      and      b'    Into  the   above 
gives   the   explicit    form 


(38)  [loj^cD^  +(co^-p^)  (coaD+  1?2^)][cd(co  cD  +  l(a^+  cd^)  )  (03(0)""'-^' 

+   leu  co^)    +   u^^^^Jl^]   +  n^u^^{(/-A)^)[{u^^-^^)^  CO 


lco2(co2-^2)^    =   0- 


This  dispersion  relation  Is  10   order  In  cu   (if 

th  2 

V   =  0,  then  It  Is  5   order  In  co  )  .   Together  with  the 

teraperature  mode,  Eq.  (33)  j  and  the  seven  modes  uncovered 

In  Section  C-IV-a,  these  frequencies  comprise  a  set  of 

eighteen  roots.   Since  the  original  relation  was  seen  to  be 

1^   order  and  all  of  the  roots  fall  into  either  one  of  the 

two  categories  considered,  it  follows  that  four  of  them  are 

covered  twice.   These  four  roots  will  be  exhibited  in  the 
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course  of  examining  some  fundamental  limiting  cases. 

(l)   Rare -Plasma  Dynamics: 

This  situation  is  obtained  by  letting  cd^  =  e^n^/e^^ 
tend  to  zero,  to  obtain 


(39)      {o?-   c^k^)^co  ^[(o)  (i  +  iTl^f  +   a^oo^]  =  0. 


The  first  factors  are  the  vacuum  electrodynamlc  roots 
(VE  in  Fig.  l).   The  second  factor  is  degenerate,  while  the 
third  is  a  pure  decay,   Ico  =  v.,  .   The  last  factor  yields  the 
four  decaying  Larmor  modes  (LR  in  Fig.  l): 


1C0+'''  =  v^  +  Ifi  t  /(n  +  Iv^)^  +  4?2^ 
(40) 

ico+"  =  v^  -  1^  t  ^(-17  +  Iv^)^  +  Kfl^ 

If  we  Impose  the  additional  constraint,   v-j^  =  0,  Eq.  (39) 
goes  to 


(i+1)     (o)^-  o^^^fo?{o^^-  n^){(/-  a^)  =  0 


where  we  have  transformed  back  to  Q,D,     from  Q,Q.,      through 

Q  =  Q  -  n 
fi^  =  fi  n 
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so  that 


(41")     fi2^2  -  (a^-  o?f=    {Q^-  o?){o?-   n^) 


(2 )   Isotropic  Plasma  Dynamics. 

The  limit  we  now  consider  Is  that  of  a  vanishing  mag- 
netic field:  fl,    n  — >  0,   wherein  we  expect  to  find  the 
classical  Interaction  of  the  plasma  with  the  VE  mode.  I.e., 
-co  +  c  k  +  03  =  0.   In  this  limit,  equation  (38)  reduces 
to 

(42)       CD^[la3  co^  +  aJ(co^-  c^k^)]^(a)  oo  +  Icu^)  =  0. 

The  middle  factor  yields  the  dispersion  relation  for 
the  phase  velocity  v  =  o/k   (MI,  In  Fig.  1.   See  below. 
Section  4-b,  for  further  discussion). 


(43)      v^  =  =_^ 


a)(a3-lv-,  j 


It  may  also  be  written  as  a  dlsslpatlve  perturbation  on  the 
more  familiar  coupling  mode,  viz., 

fhh\                  2     2     2,  2  ,  .    rCD^-c  k  T 
(44)      CD  =  CD  +  c  k  +  iv^  [ ]  . 

The  explicit  form  of  the  roots  of  the  related  cubic  are 


unreveallng,  however,  it  is  significant  that  one  of  them 
becomes  degenerate  in  the  limit  of  no  dissipation.  This 
remaining  root  behaves  as: 


(44.1)    loo  ~  v^c^lS/{o^l  +   c^k^) 


for   V-,   small.   More  precisely,  since  the  original  factor 

2 
is  squared,  two  non-zero  roots  degenerate  to  co  =0. 

The  last  factor  is  a  decaying  plasma  oscillation  with 

the  corresponding  roots. 


(45)       2co  =  -iv^  +  203^  )/l  -  (v^/2a3^)2    , 

so  that  the  purely  oscillatory  part  oo   is  decreased  by  the 
factor  /l  -  (v/2a3^)^   . 

If  the  con 
(42)  reduces  to 


If  the  constraint   v.,  =  0   is  Imposed,  the  relation 


(4  2')     ^^(^o  ■*"  ^^^^  "  '^^^^("^o"  '^^^  "  °  • 


The  last  factor  represents  pure  plasma  oscillation, 

while  the  second  represents  the  interaction  spoken  of  above. 

It  is  interesting  to  follow  these  modes  as  we  pass 

from  one  limit  to  the  next.   For  Instance,  comparison  of 

2 
Eqs.  (41 )  and  (4Z')  Indicates  that  the  two  co  =  0  roots 

2     2 
of  Eq.  (41)  have  gone  to  o^  =  cd  ,   while  the  magnetic  modes 
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became  the  degenerate  co  =  0  modes  of  equation  (45')'   Lastly 

p      O  Q  p 

the  pure  VE  mode  (cd  -  k  c  )   =0  became  coupled  with  the 

2     2  2     2  2 
plasma  in  the  above  limit  and  went  to  (cd  +  k  c   -  go  )   =0. 


(3)   Pressureless  Plasma  Dynamics. 

2 
We  now  proceed  to  examine  the  limiting  situation  cjo  — >  oo. 

The  motivation  of  this  investigation  arises  from  our 
starting  equations  (29)  and  (30).   Perhaps  the  most  drastic 
difference  between  the  equations  of  MHD  (Eq.  25)  and  these 
of  PEM  lies  in  the  Ohm's  law  equations.   While  in  the  former 
case  (viz.,   E  +  UXB  =  O),  it  serves  to  relate   E  to  U 
immediately,  in  the  latter  case  it  is  seen  to  be  a  full 
equation  of  motion  for  the  current.   Nevertheless,  comparison 
of  these  two  equations  reveals  that  they  coalesce  in  the  limit 
cOq  ->  a,.   The  momentum  equation  of  this  PEM  limit,  Eq.  (29), 
is  equivalent  to  the  momentum  equation  of  MHD,  Eq.  (l9 ) 
with  the  pressure  term  removed.   It  follows  that  in  the  con- 
sidered limit  03  — >  CO,  the  above  relation  (38)  should  reduce 

to  the  form  obtained  from  the  MHD  relation  (25)  in  the  limit 

2 
a  =0,  viz . , 
o 


(25)      v^(v^  -  a'2cos^0)(v^  -  a'^)  -  0  . 


2 
More  precisely,  we  note  that  the  Alfven  speed  A   remains  in 

2 
the  formulation  in  the  limit  co  — >  m.   It  will  be  recalled 

o 

2  2     2/ 

that  these  parameters  are  such  that  cd  ~  n  and  A  ~  B  /n 


so  that  in  order  to  maintain  a  finite  P?      as  a?      becomes 

o 

large,  the  Impressed  magnetic  field  must  also  become  large. 
The  limit  we  consider  Is  co^  -^  «>  along  the  path  A^  =  constant. 
This  constraint  Is  manifested  In  the  constancy  of  the  para- 
meter,  p   =  A  /c   =  n/co^  ,   which  Is  now  Introduced.   To 
uncover  the  nature  of  the  roots  about  large  plasma  frequency, 
we  rewrite  equation  (38)  as, 

(46)  [03^(1  +  p2)  _  p2_^2^^_^2p2  _  ^2^^  ^   p2)^ 

-^^(1  +  P)][a32(l  +  p2)  _  p2^2j  _  ^2(^2_^2)(^2_^2)| 
-  s^oi^[{o?-^^)[o?{l   +   p2^  _  p2^2^  ^  ^^2_^2^^^2^2+p2) 

_2 
where  we  have  written   e   for  03   ,  and  have  recalled 

(47)  n  =  -^  • 


Since  the  related  polynomial  Is  a  function  of  03  ,   we 

2 
will  seek  solutions  for  x  =  03  ,  of  the  form 


Note  that  the  dissipation   v^   Is  set  equal  to  zero; 
however.  Its  Inclusion  would  not  effect  the  lowest  order 
results  herein  presented. 
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(48)      X  =  x^  +  ex-^  +  .  .  .  . 

The  remaining  solutions  are  of  the  form 
(^9)     y  =  y^e"^  +  y-^  +  . . .  . 

If  these  series  are  substituted  Into  equation  (   ), 
the  series  In  x  yields  first 


(50)      x^  =  ^J^     , 

1  +  B 


which  corresponds  to  the  "cold"  fast  wave  (c"^  finite) 


fr-^s                   2    Ac 
(51)       V   -  -^ 2 


where  we  have  again  written  v   for  a/k. 


(52)       X   -   ""^ 


°    1  -H  p2   ^ 

or  equlvalently, 

2  2 
(53)      v^  =  ^  ^  ^  cos^e  =  a'^cos^0   . 
0^^+  A"^ 

Since  these  latter  four  roots  (51 )  and  (53)  are  Iden- 
tical to  those  extracted  from  the  relation  (24),  It  follows 
that  they  are  the  aforementioned  roots  which  are  covered 
twice . 
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The  series  (4-9)  yields,  to  lowest  order,  (writing  y 
for  y^) 


(5^)      ?-   y^(3  +  P^(2  +  y))  +  y((3  +  p2)(l  +  p2)  +  yP^) 


(1  +  P^)^  =  0. 


Observing  that  the  sum  of  the  coefficients  vanishes.  It 
follows  that  the  related  three  roots  are  the  plasma  frequency 
y  =  1,  together  with 


(55)  2y2^^  =  2^\   +  YP^  +  p/i+pfT+'rV'  ;   p^  -:1  +  P^ 

In  terms  of  the  relevant  physical  parameters,  these  roots 
appear  as 

(56)  20)5  ==  f^'2  +  2gd^  +  n  ivi*^   +  40)^    , 


where 


(56.1) 


and 


2t.2 
^  ^o 


(^)V(m2+  m^) 


.2    ^'^o 
(56.2)    fi*2  =^0     . 

[mVM] 
These  two  roots  (Eq.  56)  are  abbreviated  PMO '  In  Figure  1 
In  the  limit  B  — *  0,   these  roots  collapse  to 


2    2        2 
■^0  ' 


(56.3)     CO  =  o)^  ;      cd""  =  0  . 
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If  we  impose  the  constraint  Q,   »  a>      (we  have  already 
set  CD  »  all  other  frequencies),  then  Eq.  (56)  reduces  to 


(56.4)    cof  =  n^ 


the  pure  Larmor  modes. 

2  2 

To  recapitulate,  in  the  limit  co^  — >  00  ,   a  =  constant. 

the  secular  equation  (38)  tends  to 


(57)       (co^-  0)2)  (cu^- 0.2)  (0)2-  o)|)(a)2-  o)2)(oo2-  o)|)  =  0    , 


where 


(58) 


Z         p2   2  2    pVz 

O).   =  ^   ^   2    ,      03   =  2 


2     2         2       2         2       2 
03-5  =  0)^  ,      "^4  ^  ^2^0  '      ^5  ^  y^'^o 


Having  determined  the  x   coefficient  in  the  series 
(48),  the  next  order  term  x-,   is  evaluated  by  substituting 
one  of  the   x   values,  say  first  the  Alfven  mode  (52)  into 
equation  (47)  and  again  setting   e  =  0.   The  lowest  order 
terms  now  yield 

(59)       X  =  ^^    /_  p2  ^  &^  (1  +  p2^cos20 


+  ^  Ctn20(p2cos2e  -  p2) 
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which  corresponds  to  the  dispersion  relation. 


(59'  )     v^  =  -|^  cos20[i  -44  |P?  -  4  (1+  P?)cos^ 


^  ctn^eO^cos^p  -  P?)  J  ^  • 


term 

now 

yie: 

Ids 

(61) 

^1 

^0 

(-0- 

•/ 

)[1 

+ 

T 

1 

For     c   »   A,      this   reduces   to 


(60)  v^  ~   A^cos^e  \l    -  ^^   (1  +  T  ctn^e) 

CO 

^  o 

Secondly,  we  consider  the  extension  of  the  fast  wave 
roots  (50).   Again,  if  the  series  (48)  is  substituted  into 
equation  (47)  with  x^  =  ^"^J)^/!  +   p  ,   the  lowest  order 


/p^sin^e  -  cos^Qn J 
sin^e 


which  corresponds  to  the  dispersion  relation 

/^oN       2    A^c^  /,    k^c^  fi    ^  T-   /P^sin^e  -  cos^0^^ 

(62)      V   =  -^ P  ^  1 ]T-^   (1  +  -^  ( 2 >^ 

A^+   c^  [     pjo.;      ^l  sin^e 

2 
which,  in  the  limit   p   «  1,  reduces  to 


(63)      v2  ~  a2(1  -  ^  (1  -  Y  ctn^e))  . 
o 

The  singularities  about   0=0  in  equations  (59-63) 


102 


are  ficticious  and  are  due  to  the  fact  that  when  e  =  0, 
X  =  k^A^   Is  a  double  root  of  equation  (^7)  at   0=0.   It 

follows  that  the  above  forms  are  not  appropriate  about  these 

1/2 
values.   An  expression  (for  P  «  l),  good  to  0(e  ^  ),   which 

Is  valid  In  the  whole  domain  Is  obtained  by  completing  the 

square  of  the  lowest  order  term  In  (46)  and  then  changing 

variables  from  w   to  y  =  w-k;A(|i.+  l)/2,   where 

M-  =  cos  9,      and  w  =  00  /k  A  .   Retaining  terms  of  order  (  y) 

and  solving  for  w   gives 


(64)      w,  ^  4^t  V(^)^  +  |'(^^+1)(M-Y-1)  +#'(^^T-1) 


where 


Y  +  1 


e'  =  k^c%^ 


The  (+)  root  Is  the  extension  of  the  Alfven  waves, 
while  the  (-)  root  Is  the  extension  of  the  cold  fast  wave; 
In  the  direction   0  =  0,  Eq.  (64)  Indicates  that  the  two 
speeds  are  displaced  by 

(64.1)    v,„,  -v„J     =2r^=2i5Ll%m 


'"""^  °"'  '9=0  ""o  ,  m  m 

If  m  =  m,   then  the  two  waves  again  coalesce  at   0  =  o. 
The  forms  (60)  and  (63)  are  a  more  accurate  description 
of  waves  In  the  direction   0  ~  it/2.      In  this  direction  the 
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cold  fast  wave  is  decreased  by  the  factor  (l  -  k^c^/2cD^) 


(4)   Directional  Wavt 


(4-a):   k^  =  k^  ,    (k  X  B^  =  0) 


Substituting  these  values  into  the  form  (35)  yields 


(65)       F^F_(icD^  +  CO  co)  =  0: 


where 

(66)       F^  =  (k^c^-  a)^)[a)  ^  Q)  (co  t  ^)    +   iv^co]  +  co^co^ 


and  we   have   recalled  that 


(66.1)  00^  t  i^oj  -  ?^^   =    (co  jjl  n)(cD  t  ^) 


(Note  that  these  D's  are  all  positive.) 

The  root   F^  =  0,   yields  the  dispersion  relation 

(67)       v^  =  ^ . 


1  - 


(co  +  Q)  (co  t  ^)  -   iv-,co 

In  the  limit  of  m  »  m  and   v-^  =  0,   this  reduces 
to  the  more  familiar  relation,^ 

2         c^ 
(68)      v^  =  ^ 

^  "  (03  +  n) 
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2 

The  remaining  root  co  cd  +  Icd  =0,   is  the  same  root 

as  was  found  in  the  vanishing  magnetic  field  limit  (case  (2): 
Pure  Plasma  Dynamics),  and  is  the  last  factor  of  equation  (42), 

(4-b):   k^  -  0  ,    (k-B^  =  0). 
In  this  limit  the  form  (35)  reduces  to 

(69)         GD[6D(a3^-  k^C^)  +  103  03^]  ([(cD^-  k^C^  )  (o)  ^    +     ifi^  ) 

+  iaj^cu^][a)  cB  +  ifi^  +  icu^]  +  a)^n^(a3^-  k^c^)|  =  0  . 

The  first  two  factors  are  independent  of  the  magnetic  field, 
and  are  associated  with  currents  flowing  parallel  to  the 
magnetic  field,  (indeed,  they  have  already  been  encountered 
in  the  B  =  0  limit,  see  above  Eq.  42).   Specifically,  the 

non-degenerate  second  factor  is  a  generalization  of  the  so 

25 
called  magneto-ionic  root  of  the  ionosphere  theory. 

The  remaining  sixth  order  relation  represents  the 

coupling  of  the  plasma  to  the  magnetic  field,  and  yields  the 

following  dispersion  relation  for  the  phase  velocity  (o/k) : 


(70) 


where 


2 
c 


.    2      - 

lOD    CD    03 
O 


(71)  W     =     (03    03    +     i?^^)(03    03    +     iTl^    +     103^)     +    oi^    9^ 
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In  the  limit  of   v-,  =  0, 


(72)      w  =  (cD^-  a:^)(cD^-  co^) 


(73)      2a:^  =  20^  +  V?  +  o?^  t   /n^(fi^+  4fi^ )  +  cD^(2f2^  +  o)^) 


roots  Indicates  a  more  Introcate  coupling  between  plasma 
and  magnetic  field  t 
which  Is  as  expected 
In  the  limit  of 
the  pure  Larmor  modes. 


2    2 
and  magnetic  field  than  we  observed  In  the  k  ==  k  case 


In  the  limit  of  cd  — >  0,  the  above  four  roots  go  Into 


(5)   Non-Propagating  Waves. 

In  this  final  limit  we  consider  -p  — ^  0  along  with  k, 
so  that  the  resultant  phase  speeds  far  exceed   c.   The 
determinant  (35)  reduces  to 


(74)      0)^(03  CD  +  lco^)[  (03  o)  +  Ifi^  +  103^)^  +  n^co^]  =  0 


The  first  two  factors  were  already  encountered  in  the 

first  class  of  vacuum  electrodynamics.   The  roots  of  the 

(23) 
remaining  quartic  appear  as,^ 


(75)      2cD. 


^+  =  -iv-^  -  fi  t  /(Iv^  +  Vif   +  4(0)^  +  ^2) 
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'^  +  n  +   /(iv^  -  af  +  ^{o^l  +  n^) 


where  we  recall  that  the  augmented  Larmor  frequencies  Q, 
and  fi  are  related  to  the  pure  Larmor  frequencies  through 

n  =  n  -  n 

(41-) 

n^  =  n  Q,    . 

All  four  of  these  modes  are  decaying,  with  a  disslpa- 

tlve  relaxation  time,  which  for   v^  <<  cd^  ,      Is   (v-,/2)~  = 

1      o'       ^    y     ' 

The  oscillatory  part  of  the  frequency  In  this  approximation  1^ 
given  by  Re((Jo)  , 


Re(cD+)  --fit  jff2"V^ 
(76) 

Re 


(cD  )  =  +  n  t  }^ 
where  we  recall  that   H   includes  the  reduced  mass  7^   , 
(56.2)    fi*  =  eB/;^  . 

and  similarly  for  the  plasma  frequency  <x>   , 
(78)      -^  =  eV-o^  • 
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D.   Conclusions. 

I.   General  Statements. 

We  have  seen  the  dlsslpatlve  plasma  oscillatory  phenomena 
emerge  from  a  self-contained  microscopic  analysis.   This 
analysis  Indicated  that  there  are  two  fundamental  regions  of 

excitation.   These  regions  are  governed  by  two  distinct  types 

(l)  2  f2 ) 

of  dispersion  relations:   oj  =  kco^  ^  +  k  co^  '   +    ...,    and 

CO  =   cd(°)  +  kcD^^)  +  k^cD^^''  +  ...,   0)^°'*  ^   0,    respectively. 

The  first  region  Is  Indicative  of  finite  wave  propagation, 

(  11m  (o/k)  =  constant),   and  yields  Lundqulst's  Magneto 
k^-o 
Hydro-dynamics  (MHD),  wherein  the  plasma  exhibits  the  fluid 

dynamical  property  of  sound  propagation.   In  the  second  region, 

11m  (o/k)   Is  unbounded.   In  this  domain  of  phenomena,  the 
k— >  o 
plasma  acts  cooperatively  and  In  an  ordered  manner  with  respect 

to  the  external  B   field,  In  addition  to  exhibiting  an  Inter- 
action of  these  oscillations  with  the  electromagnetic  fields. 
These  Plasma  Electro-Magnetic  (PEM)  phenomena  are  also 
characterized  by  fluctuations  whose  time  variations  far  exceed 
their  spatial  variations.   The  wave  number  k  was  kept  In  the 
lowest  order  (viz.,  k  ^  O)  terms  by  Introducing  the  parameter 
P  =   ck.   If  this  parameter  Is  kept  finite  while  all  other 
space  derivatives  are  set  equal  to  zero,  then  the  analysis 
Includes  time  variations  which  are  large  compared  to  only  the 
space  variations  originating  In  the  particle  dynamics  part  of 
the  system  (i.e.,  Boltzmann  equations).   In  the  MHD  region. 
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space  and  time  derivatives  are  seen  to  be  of  the  same  order. 

The  lowest  order  (in  k)  form  of  the  master  determinant 
was  seen  to  be  1^   order  In  co.   In  the  above  scheme  of 
ordering  these  roots,  seven  emerge  In  the  MHD  limit,  while 
eleven  emerged  In  the  PEM  limit.   The  four  roots  which  were 
covered  twice  were  discovered  to  be  the  cold  fast  waves  and 
the  Alfven  waves. 

In  terms  of  the  equations  governing  the  two  regions.  In 
the  momentum  equation  the  distinction  resides  in  the  gradient 
of  pressure  term,  which  vanishes  In  the  PEM  limit  but  is 
finite  in  the  MHD*~  limit.   A  sufficient  condition  that 
Induces  such  a  separation  is  that  the  equilibrium  tempera- 
ture  T   be  either  zero  or  finite  respectively.   The  distinC' 
tion  between  the  two  current  equations  is  that  in  the  MHD 
case  it  is  the  simple  vector  relation  E  +  UXB  =  0,   while 
in  the  PEM  region  it  is  a  full  equation  of  motion  for  the 
current.   Having  observed  that  in  the  latter  equation  the 
plasma  frequency  multiplies  the   E  +  UXB   term,  it  follows 
that  the  PEM  system,  collapses  to  the  pressureless  MHD*'   as 


00  — >  (X)   (together  with  the  constraint   B  — ^  c»).   This  was 
seen  to  be  so,  with  the  modified  Alfven  speed  A   =  c  a'/(c  +  A^ 
appearing  in  place  of  the  standard  A  =  '^J^    \       •      '^^^  Alfven 

t   These  statements  are  made  for  the  related  first  order 
(in  k)  theory. 

T  The  star  superscript  on  the  MHD  abbreviation  indicates 
that  the  displacement  current  Is  included. 
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and  fast  waves  were  continued  to  terms  of  0(a)~  ),   and  a 
velocity  spread  In  the  parallel  direction  (i.e.,  to  B  ) 
between  the  Alfven  and  fast  waves  was  found.   This  Increment 
vanishes  In  the  limit  of  equal  masses.   In  the  normal  direc- 
tion the  cold  fast  wave  was  found  to  suffer  a  loss  In 
velocity  (cf.  Eq.  C-IV:  65). 

The  fact  that  these  waves  are  recaptured  In  the  limit 
k  =  0,  -^  /^  0,   Is  Indicative  of  their  electromagnetic  origin. 

In  the  PEM  limit  the  energy  equations   yield  a  decaying 
finite  temperature  difference  to  the  equilibrium  value  AT  =  0, 
with  a  related  decay  time   (2Vp)~  .   However,  in  the  MHD  class 

AT   is  zero  to  terms  of  order  K.   The  addition  of  the  MHD 

2    10     .   ** 
energy  equations  yield  a  sound  speed  a  =  -^  kT  /M. 

The  PEM  roots  in  the  limit  of  vanishing  number  density 
went  to  one  degenerate  cu  =  0,   a  purely  decaying  mode 
10)=  V-,  ;   the  vacuum  electrodynamlc  roots  (00  -  p  )   =0,   and 
four  decaying  Larmor  roots  (cf.  Eq.  C-IV:  39)'   The  fact  that 
the  Larmor  roots  remain  in  the  limit  of  vanishing  number 
density,  while  the  plasma  oscillations  collapse  to  zero  is 
indicative  of  the  fact  that  the  latter  effect  is  a  coopera- 
tive one,  while  the  Larmor  resonances  are,  as  already  stated, 
an  ordered  response  to  an  external  Influence. 


The  terminology  is  more  convenient  than  accurate.   A  more 
precise  name  would  be  "non-adlabatic  equations". 

Note:   In  limit   m  =  m,   M  =  2m  and  \^~^%   kT^m, 
cf.  footnote,  p. 
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In  the  opposite  extreme  of  vanishing  B^   field  the 
plasma  acts  collectively,  yielding  decaying  plasma  oscilla- 
tions with  a  decay  frequency,   Im  2a3  =  -v-^  ±   Im.y^cD^  -  v-j^ 
and  also  exhibits  an  Interaction  of  this  collective  behavior 
with  the  electromagnetic  fields  to  yield  six  decaying  roots 
(cf.  Eq.  C-IV:  42).   If  the  dissipation  is  set  equal  to 
zero,  then  two  of  these  roots  degenerate  to  co  =  0,  while 

the  remaining  four  to  to  the  classical  magneto-Ionic  rela- 

p    p     2  2  2 
tlon  (gl)  -  CD  -  k  c  )   =0.   The  decaying  plasma  oscillations 

become  purely  oscillatory. 

As  expected,  the  decaying  plasma  oscillations  are  also 

present  in  waves  which  are  directed  parallel  to  the  magnetic 

field.   In  addition  to  this  second  degree  root,  the  rem.ainlng 

eighth  order  root  in  this  direction  splits  into  two  fourth 

order  roots,  one  being  obralned  from  the  other  through  the 

reflexion:  +0.  -^  -Q,      -Q,  -^ +0..      The  form  of  these  roots, 

say  P^  =  {(k^c^  -  co^)  [  (cu-^)  (oM-fi)  ]  +  o^^^^  |  ,  (v-j_  =  O), 

indicates  how  the  Larmor  resonance  decouples  the  plasma  oscll- 

2     2  2 
lations  from  the  vacuum  roots  co  -  k  c  .   If  this  root  is 

Inverted  and  solved  for  v  =  cd  /k  ,   one  sees  the  genesis  of 

the  magnetic  part  of  the  MHD  roots;  for  in  the  limit  of  n^, 

~2     2  2  2      !  2 

B  — >  00  (again  writing  Q.     =  cd  6  ),  we  see  that   v  — ^  A '  . 

The  other  F_  factor  gives  the  other  v   =  A'  ,  so  that  both  the 

fast  waves  and  the  Alfven  waves  are  recaptured  in  the  parallel 

direction. 

The  decaying  magneto-ionic  roots  of  the  vanishing  magnetic 
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field  limit  are  also  present  In  waves  directed  perpendicular 

to  the  magnetic  field  (k«B  =  O).   The  remaining  sixth  order 

2     i2 
root  In  this  direction  yields  only  the  fast  wave   v   =  A 

In  the  MHD  limit.   This  Is,  of  course,  as  expected  since  the 

Alfven  waves  vanish  normal  to  the  magnetic  field. 

The  totally  non-propagating  limit   k  =  0,  -^  =   0 ,    yields 
purely  oscillatory  modes,  the  most  Interesting  of  which  are 
four  mixed  plasma-magnetic  fluctuations  (cf.  Eq.  C-IV:  75)- 
Besides  a  degenerate   oo  =  0,   the  only  other  mode  contains 
the  two  decaying  plasma  oscillations.   It  should  be  noted 
that  these  roots  are  the  formal  lowest  order  roots  In  the 
expansion  co  =  o^^°'   +   koo''  '   +  ...,   with  all   k  (including 
those  originating  In  Maxwell's  equations)  explicit.   Although 
such  roots  are  associated  with  a  rigidly  homogeneous  medium, 
the  non-lsotroplc  structure  of  the  system  Is  present.   It 
manifests  Itself  In  the  related  system  eigenvector  X^   (cf. 
Eq.  C-I:  22),  the  pertinent  components  of  which  are  (J^^,  Jj_). 
It  Is  readily  established  that  the  plasma  oscillations  are 
associated  with  the  vector  (j^^,  O)  while  the  remaining  mixed 
plasma-magnetic  oscillations  are  eigenvalues  with  the  related 
eigenvector  (O,  Jj_). 

More  generally,  as  mentioned  In  the  Introduction,  In  any 
of  the  above  limits  wherein  B^  ^  0,  those  modes  which  do  not 
contain  the  magnetic  field  are  the  eigenvalues  of  the  eigen- 
vector (J  ,  0),  while  the  complementary  non-degenerate  roots 
are  the  eigenvalues  of  (O,  J ^) . 
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The  non-lsotropic  nature  of  the  system  becomes  evident 
at  an  early  stage  of  the  analysis,  and  Is  manifested  In  the 
distinctly  different  nature  of  the  perpendicular  and  parallel 
components  of  the  wave  number,   kj_  and  k  ,  respectively. 
This  non-lsotroplc  quality  of  the  problem  permitted  the 
Larmor  resonance  phenomena  to  be  examined  In  the  limit  of 
k^^  =  0,  while  terms  of  o(k^  ~  )  were  retained  (N  =  order  of 
resonance).   This  class  of  roots  Is  four-fold  Infinitely 
denumerable,   and  Is  seen  to  Involve  0(kj_  ~  )  variations 
above  the  classical  Integral  multiples  of  the  Larmor  frequencies, 
The  first  order  (in  k)  decay  time  is  as  the  inverse  self- 
collision  frequency,  and  is  Independent  of  the  resonance  order 
of  the  mode.   The  dispersion  coefficient  Includes  the  factor 
T  (cf.  Eq.  C-II:  30).   The  form  of  this  coefficient  indicates 
that  the  dispersion  vanishes  as  B  — >  co,  while  if  n  — >  oo 
the  dispersion  is  seen  to  be  enhanced.   This   t   factor  also 
includes  the  influence  of  the  "other  gas"  on  the  dispersion 
of  the  Larmor  resonances  of  the  first  gas  (e.g.,  that  of  the 
ions  on  the  electron  resonances).   In  the  absence  of  dissipa- 
tion (v  =  0),  this  factor  (viz.,  (a(a3j^)-l)  for  the  case  of  ion 
resonance)  appears  as  N  +  m/m  which  does  not  vanish.   This 
is  indicative  of  the  difference  in  charge  of  the  two  species, 
so  that  for  any  given  N  the  charges  are  rotating  in  opposite 


More  precisely,  including  the  image  singularities  about 
a  =  -IN,  this  class  constitutes  an  eight-fold  infinity 
of  roots. 
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senses. 

The  connection  between  the  MHD  and  PEM  regions  and  sub- 
regions  Is  depicted  In  Figure  1. 

Returning  to  the  conclusions  of  the  general  analysis. 
In  the  formulation  of  the  transfer  equations,  non-local  (i.e.. 
In  time)  terms  were  obtained  which  collapsed  to  the  more  stan- 
dard local  forms  after  a  few  collision  times.   More  directly 
the  related  coefficients  are  such  that  they  become  constant 
for  values  of  the  self-collision  frequencies  which  far  exceed 
the  phenomena  frequency  oj.   In  addition,  when  the  moment 
equations  were  constructed  from  the  stated  model,  a  technique 
of  generating  higher  order  variables  was  uncovered,  according 
to  which  the  J-flux  tensor  R^  .  was  defined  (cf.  Sec.  B-II-c, 
and  App.  F).   This  latter  conclusion  was  made  evident  by  the 
lack  of  constraints  Imposed  upon  the  mass  ratio. 

In  this  vein  the  parameter  TT  Is  significant.   We  recall 
that  It  Is  the  difference  of  the  two  absolute  Larmor  frequen- 
cies so  that  If  m  =  m,  fi     vanishes.   The  difference  between 
Tl  ^   0,   /^   0     has  manifestations  first  In  the  current  equation 
wherein  we  find  that  If  Cl  =  0 ,    the  Hall  current  term 
vanishes.   Secondly,  In  the  PEM  dispersion  relation,  the  fac- 
tors which  go  over  to  the  MHD  roots  are  also  obtained  If  the 
constraint   m  =  m  Is  Included.   Physically, the  difference  in 


If  the  two  charges  were  of  the  same  sign,  then  the  disper- 
sion of  the  (rn/m)th  order  ion  resonance  would  vanish,  due 
to  the  support  of  the  first  order  electron*'resonance . 
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mass  term  Imposes  an  Inertlal  asymmetry  in  the  system.   In 

-2 
the  next  order  in   e  =  cu   ,   wherein  the  Q.     term  enters,  the 

cold,  fast  wave  is  seen  to  loose  its  Isotropic  structure  (of. 

Eq.  C-IV:  63)  and,  in  addition,  if  m  =  m  the  cold  and  Alfven 

waves  again  coalesce  in  the  parallel  direction. 

In  this  latter  part  of  the  conclusion  section  perhaps  it 
is  appropriate  to  state  some  rules  of  thumb,  so  to  speak, 
pertaining  to  the  relevant  equations  and  dispersion  relations 
of  plasmia  physics.   With  regard  to  the  governing  macroscopic 
equations  we  have  noted  that  in  the  MHD  theory,  of  the  two 
equations  --  the  momentum  and  current  equations  --  the  former 
is  the  significant  equation  of  motion,  while  in  the  PEM  domain 
the  current  equation  is  the  relevant  equation  of  motion. 
Recalling  that  these  two  equations  are  equally  significant 
since  they  are  completely  equivalent  to  the  two  individual- 
component  momentum  equations  (and,  indeed,  come  about  from  a 
change  of  variables  from  individual  flow  variables  to  total 
vluld  variables),  it  is  not  surprising  that  if  two  theories 
evolve  from  this  system,  and  in  one  of  them,  say,  the  addi- 
tion equation  (i.e.,  the  "m.oment"  equation)  is  significant, 
then  certainly  the  other  theory  should  describe  phenomena 
governed  by  the  difference  equation  (Chm's  law). 

We  have  also  noted  that  the  wave  num.ber  enters  differently 
into  the  particle  dynamics  part  of  the  system.  (Boltzmann  equa- 
tions) as  compared  to  its  roll  in  the  field  dynamics  part  of 
the  system-  (Maxwell's  equations).   In  the  latter  domain  it 
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was  natural  and  revealing  to  change  variables  from  k  to 
p  =   ck,  which,  in  addition,  insured  that   c  was  of  proper 
order  of  magnitude  (viz.,  0(e~  )).   Something  may  also  be  said 
of  the  dual  location  of  03.   For  the  most  part  we  have  seen 
that  the  dissipation  varameters   v   and   v-,   appear  in  the  form 
V  -  ico  and   v-,  -  ico.   It  is  observed  that  this  is  so  for  only 

those  CD's  which  originate  in  the  Boltzmann  part  of  the  system. 

2   P 
Factors  of  the  form  cd,  cd  -  P    ,    •••>      have  their  origins  in 

the  field  equations. 


The  asymmetry  in  the  perpendicular  and  parallel  (i.e., 
to  B  )   parts  of  the  formulation  was  Immediately  revealed  by 
the  expansions  involved  in  the  analysis.   These  were  seen  to 
divide  into  two  distinct  classes. 

Firstly,  the  expansions  about   k   were  present  in  the 

development  of  the  S  C   functions.   These  forms  are  functions 

of   k   through  the  parameter  a  =  (v-ico  +  ix)/fi,   where 

X  =  k  ^   .   The  related  expansions  involve  terms  of  the  form 
II  11 

--1   _     >^l  /-I    _        "    "    ,  \ 

^        ~    v-io)  ^     v-icjo    •  •  •  ^ ' 


Q 

The  expansion  parameter  is  ^/  _j_^^  x  =  ex,   where   x   is 
non-dimensional  through  the  wavelength  A,   and  the  thermal 
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speed  C,   and  Is  of  0(l).   In  the  MHD  limit,  v  »  co,  so 
that  £  ~  '>\  /^,      where   A   is  the  mean  free  path,  A  ~  c/v. 
In  this  limit  the  related  series  In  x   Is  seen  to  converge 
rapidly  If  A  »  A^ . 

In  the  PEM  limit,  wherein  gd  remains  finite  as  k  — >  0, 
it  is  relevant  to  consider  cd  »  v.  In  which  case  e  -  C/v, 
(where  v  =  Acj) ,  so  that  the  convergence  condition  becomes 
V  »  C,  which  in  this  limit  is  usually  the  case. 

The  second  class  of  expansions  are  centered  about  small 
kj_,  and  appear  in  the  development  of  the  exponential  form 

e       =1  +  1 


n 


The  relevant  expansion  parameter  is   k|/n  =  Ly/A ,   where   y 
is  non-dlmenslonalized  through  the  Larmor  radius   L  =  C/Q,, 
and  the  thermal  speed  C,   so  that   y   is  of  0(l).   The 
condition  for  rapid  convergence  now  appears  as   A  »  L. 

The  criterion  to  be  satisfied  in  order  that  the  above 
sufficient  conditions  (viz.,  A  »  L,  A  )  be  realized,  is  that 
these  lengths  be  small  compared  to  the  dimensions  of  the 
environment.    To  these  ends  we  will  now  compute  the  mean 
free  path  and  the  Larmor  radius  for  a  few  relevant  physical 
environments.   All  of  the  following  are  order-of -magnitude 


We  note  that  these  are  only  sufficient  conditions.   Indeed, 
many  of  the  results  which  were  herein  obtained  have  been 
observed  in  the  laboratory  where  Aq ,  in  most  instances.  Is 
gigantic  compared  to  the  apparatus  lengths. 
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_2h, 
calculations.   The  mass  Involved  is  taken  to  be  10    gms , 

(~  proton  mass) . 

The  mean  free  path  A.,   is  given,  approximately,  by 


The  distance  d  is  the  Debye  distance 


d2  =  1^  =  t^  ,  d(cm)  .  7(Vn^)l/2 

^  "o    "^o 


so  that  for  A   we  may  write 


A^  =:  3  X  10^  T^n^  (cm) 


The  Larmor  radius  is  given  by 


L(cm)  =  C/(^  ^  (V5)Ty2/B^ 


The  following  table  has  been  constructed  with  the  use 
of  these  formulas. 
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The 

Sun'"' 

0<>) 

Sherwood 

Inter- 
stellar 
Space^^^) 

Center 

Photo- 
sphere 

Chromo- 
sphere 

Corona 

Units 

^0 

108 

10^-102 

10^ 

10* 

lo*-lo5 

io*-io6 

°K 

^o 

10^5 

10° 

102* 

10^7 

10^° 

10^ 

cm-^ 

Bo 

10* 

10-5 

10 

* 
10 

* 
10 

* 

10 

gauss 

^C 

10* 

lO^l-lO^ 

10-^ 

10-6 

10-10^ 

10-10^ 

cm 

L 

10° 

10^-10^ 

10^ 

10 

10 

10-10^ 

cm 

It  follows  that  the  related  inequalities  are  more  readily 
established  in  the  astronomical  regions. 


These  values  become  very  doubtful  in  the  inner  regions 
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D-III.   Related  Problems. 


We  have  just  noted  that  the  included  formulation  was 

constructed  about  the  expansion  of  the  s'^c'^  functions  and 

the  e""""^  function.   It  will  be  noted  that  in  the  event  that 

B  becomes  small,  the  latter  expansion  Is  poorly  convergent. 

The  related  formalism  which  includes  the  small   B   domain 

o 

phenomena  is  uncovered  if  we  recall  the  form  of  the  moments 
1^,     (cf.  Eq.  C-II:  l), 

(The  notation  is  self-contained.)   We  are  now  interested  in 
the  expansion  of  this  integral  about  large  a.   Such  an  expan- 
sion is  generated  by  a  parts  operation  of  the  following  form; 


^e^/Cdf  =|)e^a'  ^  df  =  ^  [e-]'[l]d)/-  e^  §^   j  -^e^[^]' 


2lT 

0 


In  this  manner  we  obtain  the  desired  asymptotic  serle; 
(recalling  a  '■-^e~  ),  which  to  third  order  appears  as 

In  the  PEM  investigation,  this  difficulty  was  obviated  when 
Bq  — >  0,  since  the  analysis  was  such  that   A  was  a  higher 
order  infinitesimal  than  L   (i.e.,  it  vanished  first  in  the 
limiting  process). 
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Let  us  consider  the  lowest  order  term, 

f  ^.  exp|2a'i/'+ 2ik    -^  ,  sin  "^f     / 

(i)  e^Ade £ -^ ^(^) 

J  2a  +  2ik_L-|_L^°s  '^ 


27r 


Recalling  that 


is  periodic  in  "ii ,      we  obtain 


is!!L^A{o) 


2a  +  2ik^-^j_ 


The  total  moment   I 


appears  as 


(   d^^  ^   h(^)H^(^) 

V  -  ico  +  Ik"^ 


This  first  term  is  recognized  to  be  the  I        moment  of  g 

with  B  =  0. 
o 

Another  source  of  remaining  endeavor  is  found  in  the 
integrations  over  4   (cf.  p.  67).   We  recall  that  the  region 
of  investigation  in  the  cd  plane  was  restricted  to  Re(-icD)  f   v 
in  order  to  maintain  well  defined  integrations.   It  is  at 
this  point  that  a  more  intensive  transform  investigation 
would  uncover  the  influence  of  the  dissipative  parameters 
and  the  Larmor  resonances  in  the  extended  Landau  problem. 
It  is  interesting  to  observe  that  such  a  formulation  is  con- 
structed about  the  parallel  part  of  the  analysis,  which  would 
imply  that  even  these  extended  phenomena  with  an  impressed 
magnetic  field  still  have  their  origins  in  the  electric  part 
of  the  fields. 
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In  conclusion,  the  standard  problem  remains  of  the 
extension  of  the  related  expansions  about  large  wavelength 
to  higher  order  terms,  together  with  the  complementary 
Investigation  about  small  wavelength,  which  would  start 
along  the  lines  of  the  above  series  in  a 
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Appendix  A:   The  Collision  Forms. 

In  this  appendix  we  will  show  how  the  postulated 
equations  (B-II:  1,2)  evolve  as  a  first  step  In  a  certain 
Itteratlon  scheme  solution  to  the  full  Boltzmann  equations. 
The  procedure  Is  such  that  the  correction  to  the  local 
Maxwelllan   f   Is  given  by 


f(N)  ^ 

^0  + 

Si  - 

or 

equlvalently. 

f(N)  ^ 

,(N- 

1)  + 

These  correction  terms  are  determined  by  the  Boltzmann 
equations,  written  In  the  form  (say,  for  f, ): 


^^-H 


Df^  =  /  [ff -1)  fj^-1)   -  ff )  fJN)] 
^  J    [f(N-l)V(N-l)'  _  f(N)  f(N)^ 

+  —  +        —  T 

The  operator  D  represents  the  convectlve  derivative, 
D  ^  |^+  ^  V  +  (e/m)^  (E+^XB).V| 


The  notation  In  these  appendices  Is  Independent  of  the 
preceedlng  analysis  and  wherever  necessary  Is  redefined. 

It  has  been  shown  that  these  forms  also  follow  from  an 
elgenfunction  expansion  of  the  linearized  collision  Inte- 
grals (H.  Grad  and  L.  Slrovlch,  private  communication): 
cf.  also  ref.  Ol) ,    Sec.  30,  p./24-,  and  refs.  (3^,  35). 
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and 


d^ij^^J  IV,.  Ida,. 


V    is  the  relative  velocity  (^.-  ^.).   The  Coulomb  crosf 
section  Is  given  by 


Y^.   dc^.   =  [e^/kTreJ?7f{ci-Y)-^dn 


where  a   Is  a  unit  vector  that  bisects  V  and  V .   The 
prime  and  subscript   1   notation  are  standard  (cf.  references 
19^  31)-   The  element  of  solid  angle  is  about  the  spherical 
angular  coordinated  (0,  )6)    of  the  vector  a,   so  that 
(0  <  0  <  |). 

Let  us  consider  the  first  Itteratlon,  calling  f^  '  =  f, 
and  dropping  the  (+)  subscript  on  the  functions  in  the  self- 
collision  terms. 


Df 


dS 


+  =  f{^o  ^;,i  -  [ Vo,i  +  Vi  +  ^0,1  s]}^-i-+ 

-  /  (^o^l  +  ^0,1  g)^V  +  J+J^o^^^  ^o,  J 
/  (^oS--*-  ^o,-S+)^^+-  ' 


neglecting  non-linear  terms  in  g.   The  operator  J_^_  repre- 
sents the  cross  collision  integral. 

For  Maxwellian  molecules  we  recall  that  Vda   is 
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Independent  of  V   so  that   /  e.i<^_^_^.   =  /  S_d2_^_  =  0,   as 
these  terms  Include  the  factor  /  (f^-  f)d  I  =  n  -  n  =  0, 
recalling  that  the  first  three  moments  of   f^  are  equal  to 
those  of   f.   If  we  continue  this  property  to  the  case  of 
Coulomb  particles  the  above  becomes 


D^+  =     ^  t^O-  ^]  -^  J+J^O,4-'  ^O,  J  ' 


The  desired  form  follows  (i.e.,  with  the  appropriate  simpli- 
fications imposed  on  the  collision  frequencies),  if  these 
collision  terms  are  linearized  about  the  absolute  Maxwellian 
states  (f|°);  f(°)). 


Appendix  B:   The  Formal  Expansion. 

The  algebraic  system  of  transform  equations  (C-II:  l) 
are  written  symbolically 

(1)       A. .  X.  =  0. 

A   is  a   10  X  10  matrix  (with  the  thirteen  external 

parameters:   e,  e  ,  M-  ,  m,  m,  B  ,  n  ,  k,  T  ,  v,  v,  v-,  ,  v^j, 
■^  o   o    ^   "^   o   o      o         1   ^ 

and  X  is  a   10  x  1   collumn  vector  the  components  of  which 
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are  transform  perturbations:   [n,  u,  T;  n,  u,    T] . 

The  technique  of  solution  Is  to  expand  these  variables 
In  a  power  series  about  k  =  0  (where  k  Is  appropriately 
non-dlmenslonallzed) ,  so  that 

(2)  03  =  a)(°)  +  ko^^l)  +  k2a.(2)  +  ..., 
so  that, 

(3)  A(k,a))  =  A(k,  03^°^  +  ko)^^)  +  ...) 


=  A(0,  J''^)+(A^  +   A^a3^^);k 


Similarly,  for  the  eigenvector  X  we  write 


(K)  X  =  x(°)  +  kX^l)  +  k2x(2) 


If  these  latter  expansions  are  substituted  Into  AX  =  0, 
and  coefficients  of  equal  powers  In  k  are  equated,  there 
results 


a(°)(0,  a3(°))x(°)  =  0 


a(°)(0,  a)(°))X^^)+  [A^(0,  J°)-':.J^hjO,J''hu^°^    -   0 


This  system  may  be  rewrltted  In  the  form 

-  126  - 


(6)      a(°)  x(°^  =  0 

(6.1)      a(°)  X^^)  =^  J^{o:>^^h    ',  |A^°)|  =  0 


(6.2)      a(°)  x(2)  .  y2(a.(2))  ,         |a(o)|  .  q 


The  m   equation  appears  as 

(6.m)      a(°)  X^"^)  =  yM""^)    '         |A^°^I  -  0   . 

The  vanishing  of  the  determinant   Ia^'^'^I   Insures  the 
existence  of  non-trivial  solutions  to  equation  (6),  and 
serves  to  determine  0)^°'',   whence  X^°''.   This  condition, 
I A    1=0,   Implies  that  the  linear  equations  comprising  (6) 
are  not  independent  or,  equivalently,  that  the  column  vectors 
comprising  k^° ^    are  linearly  dependent.   It  follows  that 
solutions  to  (6.m)  will  exist  only  when  the  components  of  the 
vector  y   are  linearly  dependent  in  the  Identical  manner  as 
are  the  column  vectors  of  A^^-^.   Formally,  for  matrix  equa- 
tions such  as  (6.m)  this  necessary  and  sufficient  condition 
for  the  existence  of  solutions  is  that  all  the  left  null  vectors 


of  A^  '^   be  orthogonal  to  y  .   (The  more  general  theorem 


covering  the  case  when  A^  -^   is  an  arbitrary  linear  operator 

is  that  the  right  null  vectors  of  the  adjoint  of  A^  -^   be 

56 
orthogonal  to  y  • ) 
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Geome 


trlcally,  if  A^°'^   Is  of  rank  r  and  dimension  n. 


It  maps  all  vectors  In  R   onto  the  hyper space   R  ,   so  that 
for  A*'°''X  =  y  to  be  meaningful,   y  must  lie  In  R  .   The 
normal  vectors  to  this  hyperspace  are  the  left  null  vectors 
of  A^^''  (viz.,  If   z   Is  such  that  /zA^°^|  =  0,   then 
<^z|a(°)x)=<;za(°)|x)=  0). 

Returning  to  equation  (6.m),  suppose   y  Is  a  left  null 
vector  of  A^°S   then  solutions  will  exist  only  If 


(6.m')    <Yiy^(a3("^))^=.  0 


which.  In  turn,  serves  to  determine  co^  '  .      Consider  that 
A^°''   Is  of  rank  r,   dimension  n  (n  =  10   In  the  original 
system).   There  will  be  (n-r)  redundant  equations  In  the  set 
(6.m)  when  co   satisfying  (6.m')  Is  resubstltuted.   Discarding 
these  (n-r)  equations,  the  remaining  system  of  r   equations 
In  n  unknowns  yields  a  (n-r)  parameter  of  solutions  X^  '' . 
Specifically,  equation  (6.1')  appears  as, 

(1)  <rU,x(°)> 

Resubstltutlon  of  this  value  Into  equation  (6.1)  yields 
X^   ,   so  that  the  problem  is  solved  to   0(k)  terras,  and 


appears  as 


CD  =  03^°)  +  ka)^^) 


X  =  x(°)  +  kx^i)  . 
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The  higher  order  corrections  are  similarly  developed. 


Appendix  C:   Other  Methods  of  Integration. 

In  Section  C-III  (p.  67)  it  was  noted  that  the  d^^^ 
integrations  remain  bounded  if  Re  (-loo)  >  -v.   An  alternate 
procedure  may  be  Introduced  in  the  early  stages  of  the 
analysis  to  cover  the  complement  of  roots:    Re^ico)  >  v. 
To  these  ends  we  recall  equation  (C-ll:  l)  for  the  i,        moment 


of  g  and  consider  specifically  the  magnetic  weighting  factor 
(e^^^-  1)"^.  In  the  event  that  Re (a)  <  0  (which  is  satisfiec 
if  Re(ico)  >  v),   the  above  factor  may  be  expanded 


[1  -  exp(47ra)]  ^  =  ^ 
n=0 


n^Tra 
e 


so  that   I„  becomes 

I   =  -2    <5?^^        6.1/   exp|2at  +  2iy  sin  ^^ } -^H^g  . 
o 

In  this  scheme  of  Integration,  three  of  the  remaining 

four  integrations  are  readily  performed,  while  the  remaining 

integral  is  not  expressible  in  terms  of  well  known  functions 

(except  for  certain  specified  directions  of  propagation). 


The  included  analysis,  of  course,  also  covers  this  comple- 
mentary region;  cf.  footnote  p.  67. 
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Another  method  of  obtaining  a  more  Inclusive  set  of 
formulas  Is  obtained  If  we  rewrite  the  factor  e^^V(e  ^^-  l) 
In  terms  of  Its  Fourier  representation: 


e2aV(e^^^-   l)    = 


27r(2a   -   In) 
The   Integral,    C-II:    1,      now  appears   as 


d^^  (^dif  ^  exp[lm/^  +21n  sln-^l 
7r(2a  -  In) 


t>(V^^^^^-^^) 


The  Integrations  over  f   yield  linear  combinations  of 
Bessel  functions  (cf.  p.  53)- 


Appendix  D:   Tabulation  of  Integrals. 

In  the  following  tabulation,  we  recall  the  definition 

,1  ,27r 


<^>1 


/5Fd^|  =  i 

(2Tr) 


372 


o  -1  o 


■^V2  ,  ,2 


F  I  d|  dcos 


and  H^   =  (4  -  3)/2  . 


4^^>=  1.3  ...  (2n+l);        4f)>=  1-5  •••  (2n-l) 
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hf   ef>  \     ?2("-)  T„„  ,   Where 


T 

nm 


X   (1-x  )  dx  . 


/^2m  |2n^^  1.3  ...  (2(m+n)  +  l) 
\  ^    J  ^      2(2(m+n)  +  l) 


/^2m  |2n\   1.3  ...  (2(m+n)  +  l) 
\      J  /        2n+l 

<'{e-A)?)=A^<'e^?^)  =  A^5^^  =  A 
<^(e-A)(5-C)^  =((l.-h){i-C)ii^   =  A-C 

-  k^(Bxh)xB  +  2(k_LX  B)k_|_' (BKh^)  . 


This  latter  Integral  appears  in  the  coefficients  that 
multiply  the  [S^]°,  [Pc]^^\    and  [C^S^]^°^   terms  in  the 
velocity  moment  equations  (cf.  Sec.  C-III-b). 

The  following  integration  appears  in  Section  C-II: 
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L  =<(|-V)^(k,.|j2n-2\   _1     f  ^-e/2    ^.^  ^(k,.^j2n-2  ,3^ 

where  V  Is  such  that  B-V  =  0. 

Let  us  choose  a  representation  wherein  k^  =  (kj_,  0,  O) 
and  the  cylindrical  coordinate  of  the  vector  ^  is  measured 
from  kj_   in  the  plane  normal  to  B.   It  follows  that 

^   =    [£,j_(cos  a,    sin  a);    ^^J 


e  -  ii  .  e 


(k,.ej2"-2   =   kf-2   5f -2   co.2n-2  a    . 


and 


/   L, 


V^ 


2n-2 


(27r) 


372 


dl. 


de 


2lT 


dae 


■^^/2(v^cosa  +  VySlna)|f+lcos2'^-2a 


/cos  a\ 

ysincxy 


The   three    standard   integrations   that   are    included  are: 
^   cos^^ada   =  27r  b^^,      /V'"''^   exp(-y^/2)dy  =   2^nl,      and 

/^°°exp(-z^/2)dz   =    (27r)-^/^.      If     k      is   a  unit   vector   in  the 

—  CO 

k,    direction,    then. 
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L(V)  =  2^nl  k^^  ^[^211^  ^'^  ^  ^2n^  ""  ^^^  "^  ^'"^^^ 


2^    2.i^...(2n)         2n    2n^ 


Appendix  E:   A  Note  on  Determinants. 

In  the  early  stages  of  this  Investigation,  it  was 
convenient,  when  working  with  the  two  dimensional  vector  u^ 
to  employ  a  complex  representation  according  to  which 
Uj_  =  u  +  iv,   so  that,   for  Instance,   B  X  Uj_  =  iUj_.   Due 
to  the  fact  that  the  related  vectors  are  transforms,  they 
in  turn  are  already  complex,  and  an  evident  confusion  arises 
specifically,  over  the  presence  of  the  different  nature  of 

the  two  origins  of   1.   The  answer  to  this  problem  is  that 

2 
the  procedure  is  consistent,  aid  all  i's  are  such  that  1  =-1, 

That  this  Is  so  follows  from  writing  out  the  full  coordinate 
representation  of  the  vector  equations,  and  showing  that  the 
related  n  x  n  determinant  may  be  written  as   c  c,   where 
c   In  an  (n/2  x  n/2 )  determinant  and   c   is  its  complex  con- 
Jugate.   The  determinants  which  arise  are  of  the  form 


c 
-m 

f 
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where   c   Is  the   2x2 


a  +  lb 


p  +  ie 


f  +  ig 


The  submatrices  are  the  forms  which  arise  In  the  complex 
representation.   More  generally,  the  statement  Is  easily- 
extended  to  D   of  dimension   n  (n   even). 

Formula  (l)  Is  seen  to  be  a  generalization  of  Sovlllart's 
observation-'^'  that  when  p=-e,  e=w,  f  =  a,  g=-b;  then 
D  =  (a^  +  b^  +  c^  +  d'^f. 


Appendix  P:   The  Transfer  Equations  of  a  Plasma. 

Let  the  charge  and  mass  of  the  first  and  second  component; 
of  the  plasma  be  (e,m)  and  (-e,m)  respectively,  where  e,e  >  0. 
The  two  Boltzmann  equations  appear  as 

(1)  1^  +  l-V7f  +  e/m  (E+  ^XB)-\7^f  =  j(f,f)  +  J(f,f)  =  J 

(2)  II  +  ^-Vf   -  e/m  (E+^XB)- V^f  =  j(f^f)  +  3-{f,f)    =   ~J         • 

_  -      -  — - 

This  problem  Is  considered  In  detail  In  reference  (poj. 


13^ 


The   J   terms  represent  the  collision  Integrals. 
Consider  the  change  In  variables 

me  +  me 

O) 

G  =  ef  -  ef      f  =   %-=4^   . 
me  +  me 

The  Boltzmann  equations  for  F  and  G  are 


(4)       1^  +  4- VF  +  (E  +  4XB)-\7^G  =  mj  +  mj 


(5)       If  +  ^\7G  +  (E  +  4xB)-\7^(aF  +  pG)  =  ej  -  ej 


The  constants  a   and  p   are  given  by 
(6)       a  =  ee/mm  , 


mm(em  +  em) 
so  that 

a(e,m;  e,m)  =  a(e,m;  e,m) 
(7) 

P(e,m;  e,m)  =  -^{e,m;    e,m)  . 

If  e  =  e   and  m  =  m,   p   vanishes. 

In  the  following,  <^a'^  represents  the  integral  of  A 
over  velocity  space.   The  mean  mass  flow  U   is  defined 
according  to 
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(8)  pu  =(ny  , 

where  p   Is  the  mass  density, 

(9)  P  =  6)  . 

The  extraordinary  velocity  c   Is  given  by, 

(10)  c  =  ^-u  . 

The  charge  density  Q  Is 

(11)  Q  =<'g)  . 

The  conduction  current   J   Is  related  to  the  total  current 
J '   through 

(12)  J  =^Gc)  =  \G(|-u)^  =  J'-  uQ  . 

The  pressure   P.  .   and  "J-flux"   R.  .   tensors  are  given  by 

(13)  P,.  =^Pc,c.^  ,    R,.  =(go,=  .>. 

The  Internal  energy  W  and  Its  charge  counterpart   R,   are 
given  by 

(1^)       W  =  I  \Pc^)=  I  Tr  P.J  ,    R  =  |<(gc2^  =  I  Tr  R.  ^  , 

Finally,  the  heat  flux  Z  and  the  R  flux  K  are  given  by 
(15)      Z=|<(fc^c^  ,    K  =  |'(gc^c^   . 
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The  equations  governing  these  variables  are  obtained  by 

operating  on  equations  (4)  and  (5)  with   /  d  2, ,   /  d^?,c, 

1  3   2 

77  /  d  |c  .   (in  forming  these  equations  we  have  recalled 

2  he.         h^ 


that  <Fc>=  II  =  |i  =  0,   ^  =  ^=5,..) 


The  continuity  equations: 
(16) 

The  momentum  and  current  equations: 

P{^  +   u-  V)^  +   V  -Pij  -  Q(E  +  uxB)-JxB  =  0 

(17)  (|^+  \7-")J  +  {q  It  +  (-^  +  uQ).VJu  +  T7-R.J. 

-  pa(E  +  uKB)  -  p(EQ  +  JXB)  =  <(c(ej  -  ej)^ 

The  energy  and   R   equations: 

(|^+    y-u)W  +  '\7-Z+    P    :V^    -    J-[E  +   uXB]    =    0 

(18)  (|^+    V-u)R  +  V  -K  +   R    :    Vu    -   J-[(|^  -   u.V)u 

+  P(E  +  uxB)]    =  I    <(c^(ej    -   ej)^ 


R    :    r7u    =   R.  .  V-U  • 


157 


Bibliography 

S.  Lundqulst,  Ark.  Fys.  5,  297  (1952). 

K.  0.  Frledrlchs  and  H.  Kranzer,  Institute  of  Mathematical 

Sciences,  New  York  University,  NYO-6^86-VIII  (1958). 

J.  Bazer  and  0.  Plelschman,  MH-10,  NYU,  IMS  (1959)- 

A.  A.  Blank  and  H.  Grad,  NYO-6486-IV,  NYU,  IMS  (1958). 

D.  Bohm  and  D.  Pines,  Phys .  Rev.  85,  2,  338  (1952). 
I.  Bernstein,  Phys.  Rev.  IO9,  1,  10  (1958). 

E.  Gross,  Phys.  Rev.  82.  2.  232  (1952) « 

K.  N.  Stepanov,  JETP  (USSR)  3i,  1292  (1958). 

P.  Auer,  H.  Hurwitz,  Jr.,  R.  Miller,  Phys.  Fluids  3^,  6, 

501  (1958). 

W.  P.  Allis,  Conf.  on  Plasma  Oscillations,  Union  Carbide 

Corp.,  Ind.,  Ind.  (1959). 

A.  Vlasov,  J.  Phys.  (USSR),  9:25,  130  (19^5)- 

E.  G.  Harris,  Naval  Research  Lab.  Rept .  49^4  (1957). 

H,  Grad,  NY0~25^3,  NYU  (1959). 

L,  Slrovich,  NYO-9086,  NYU  (1960). 

E.  Gross,  D.  Bhatnagar,  M.  Krook,  Phys.  Rev.  91,  3,  511 

(195^). 

L.  Spltzer,  Jr.,  Physics  of  Fully  Ionized  Gases, 

Intersclence  (1956). 

H.  Grad  and  M.  H.  Rose,  NY0-6486-II,  NYU  (1956). 

I.  Kolodner,  NYO-7980,  NYU  (1957). 

S.  Chapman,  T.  G.  Cowling,  The  Mathematical  Theory  of 

Non-Uniform  Gases,  Cambridge  University  Press  (1953)- 

-  138  - 


20.  I  Bernstein  and  S.  Trehan,  MATT-42,  Princeton  Univ.  (i960) 

21.  A.  Banos,  Jr.,  Phys.  Rev.  91,    6,  1^35  (1955). 

22.  J.  H.  Plddlngton,  M.  Notices,  Roy.  Astro.  Soc,  II5, 
671  (1955). 

23.  H.  Grad,  NYO-6486-III,  NYU  (1956). 
2^.   H.  K.  Sen,  Phys.  Rev.  88,  816  (1952). 

25.  L.  Oster,  Rev.  Mod.  Phys.  32,  1,  l4l  (1952). 

26.  L.  Landau,  J.  Phys.  (USSR),  10,  25  (19^6). 

27.  K.  M.  Case,  Ann.  Phys.  7,  3^9  (1959)- 

28.  P.  Reichel,  NYO-7697,  NYU  (1958). 

29.  B.  Zumlno,  Phys.  Rev.  IO8,  III6  (1957). 

30.  A.  Sltenko,  K.  N.  Stepanov,  JETP  (USSR)  31,  642  (1956). 

31.  H.  Grad,  "Theory  of  Rarefied  Gases",  in  Rarefied  Gas 
Dynamics,  Pergamon  Press  (1960). 

32.  T.  G.  Cowling,  Magnet ohydrodynamic s ,  Interscience  (1957). 

33.  G.  P.  Kulper,  The  Sun,  University  of  Chicago  Press  (1953). 

34.  P.  Welander,  Ark.  Fys.  7,  44  (195^)- 

35.  E.  P.  Gross,  E.  A.  Jackson,  Tech.  Note  No.  8,  Off,  of  Sci. 
Res. ,  USAF   (1959). 

36.  F.  Friedman,  Principles  and  Techniques  of  Applied  Math 
John  Wiley  (1957) • 

37.  C.  Souillart,  Nouv.  Ann.  de  Math.  XIX,  320  (i860). 

38.  R.  Herdan,  B.  S.  Liley,  Rev.  Mod.  Phys.  ^K,    731  (1960). 

39.  E.  Astrom,  Arkiv  Fysik  2,    42,  448  (1950). 


139 


NEW  YORK  UNIVERSITY 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

4  Washington  Place,  New  York  3,  N.  Y 


JAN2  6,962„^^^^,^ 

O..UO.O 

"— ~— 

/ 

r^  .     1 

NYU                                          c.l 
1      NYO-9762 

Liboff,   Richard,    L 

AUTHOR 

TIYXJ                                           c.    1           ^• 
■       ^TYO-9762 

Liboff,   Richard,    L 

AUTHOR 

Propagation  of  long  wave- 

f;TLEl<^ii8th  disturbances   in 
plasma. 

i     °^7 

BORROWERS  NAME 

NUMBER 

1] 

J 

N.  Y.  U.  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 

